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ABSTRACT 


~p  This  report  sunmarizes  progress^ during  the~~fifth~yegr~u^?fundamental 
research  on  welding  processes.  Studies  include; metal  vapors  in  the  arc; 
development  of  a  high  speed  infrared  temperature  monitor;  digital  signal 
analysis  as  a  weld  process  monitor;  convection  in  arc  veld  pools;  droplet 
transfer  and  contact  tip  wear  in  gas  metal  arc  welding  of  titanium;  and 


fracture  toughness  of  HY-80  weldments.  l\e-< 


s.-:  s  ; 


V') 


I.  INTRODUCTION 


This  report  describes  work  performed  in  the  MIT  Welding  Laboratory  under 
Office  of  Naval  Research  sponsorship.  The  work  is  generally  fundamental  in 
nature,  but  attempts  have  been  made  to  interface  the  studies  with  specific 
U.S.  Navy  Programs.  The  best  examples  of  the  directed  programs  include  the 
evaluation  of  the  fracture  toughness  of  HY-80  weld  metal  and  the  study  of 
contact  tip  wear  in  titanium  welding. 


The  research  conducted  from  15  February  1983  to  15  February  1984  can 

broadly  be  divided  into  four  main  topics,  viz.; 

V 


1. ^  Spectroscopy^nf  arc  welding  • 

2.  Welding  of -heavy  section  titanium  •  - 

3.  "^Fracture  toughness.  o^MTI-80  weld  metal 


4.  Automation  of>arc  welding  o  u  f~p /?>  &  j t  ^  f. 


II.  SPECTROSCOPY  OF  ARC  WELDING 

Spectroscopic  studies  of  the  arc  welding  process  performed  during 
the  past  contract  year  include  development  of  a  Multichannel  Infrared-Red- 
Temperature  Micro-Analyzer,  and  a  study  of  light  emissions  from  the  arc. 

These  two  studies  are  coupled  through  a  thermodynamic  upper-bound  analysis 
which  has  been  published  previously  (1) ,  and  was  expanded  in  a  recently 
published  paper  (Appendix  A) . 

A.  Metal  Vapors  in  the  Arc 

Metal  vapors  in  gas  tungsten  welding  arcs  are  being  studied  in  order 
to  determine  the  effects  of  these  vapors  on  arc  properties,  and  subsequently, 
on  weld  bead  configuration.  Emission  spectroscopy  and  monochromatic 
photography  were  used  to  determine  the  composition  and  distribution  of 
vapors  in  the  arc. 

This  project  has  produced  one  thesis,  the  abstract  for  which  is  included 
in  Appendix  B.  Two  papers  are  now  being  prepared  for  publication  based  on 
this  work.  Some  results  for  this  project  were  also  presented  at  an  AIME 
conference.  The  abstract  for  this  presentation  is  included  in  Appendix  C. 
Another  sponsor  (U.S.  D.O.E.)  has  funded  a  continuation  of  this  work  into 
an  examination  of  light  emissions  from  the  arc  as  a  weld  process  sensor. 

B.  Temperature  Measurements  of  the  Weld  Pool 

In  order  to  experimentally  confirm  theories  on  weld  pool  convection  and 
vaporization  cooling,  this  project  is  an  attempt  to  produce  absolute 
temperature  distribution  maps  of  the  weld  pool  surface.  A  unique  radiation 
pyrometric  technique  is  being  developed  to  make  these  measurements  using  a 
Multichannel  Infrared-Red  Temperature  Micro-Analyzer  (MIRTMA)  prototype 


device. 
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This  work  has  produced  one  thesis,  the  abstract  for  which  is  included 
in  Appendix  C,  and  one  technical  paper  (Appendix  DOE)  presented  at  a  SPEI 
conference.  Another  paper  has  been  accepted  for  publication  and  one  or  two 
more  are  now  being  prepared. 

Research  in  this  area  is  now  continuing  in  order  to  refine  the  temperature 
measurement  technique  and  develop  an  improved  MIRIMA  device.  It  is  predicted 
that  this  measurement  technique  will  have  applications  beyond  the  welding 
field. 


III.  WELDING  OF  HEAVY  SECTION  TITANIUM 

Joining  of  heavy  section  titanium  is  of  particular  concern  to  the  Navy 
for  advanced  systems. 

A.  Droplet  Transfer  in  Gas  Metal  Arc  Welding 

This  project  has  been  postponed  for  this  contract  period  to  fund  other 
projects  and  until  an  interested  student  could  be  found.  It  is  expected 
that  this  research  will  be  restarted  during  the  next  contract  year.  During 
the  past  year,  one  paper  has  been  published  (Appendix  E) . 

B.  Contact  Tip  Wear  in  Gas  Metal  Arc  Welding 

Copper  contact  tips  wear  very  rapidly  when  titanium  electrodes  are  used 
in  gas  metal  arc  welding.  The  goal  of  this  work  was  to  identify  the  mechanisms 
causing  the  increased  wear  and  to  find  methods  to  reduce  the  wear.  During 
the  past  year,  this  project  has  produced  one  thesis  and  was  presented  at  an 
AIME  conference,  the  abstracts  for  which  are  in  Appendices  F  and  G,  respectively 


This  project  has  been  discontinued. 


IV.  FRACTURE  TOUGHNESS  OF  HY-80  WELDMENTS 


HY-80  is  an  important  alloy  for  a  variety  of  Navy  systems.  The 
fracture  toughness  of  thick  plate  submerged  arc  welds  is  of  particular 
interest.  This  project  is  an  investigation  into  the  metallurgical  factors 
influencing  charpy  energy  of  these  weldments. 

Work  on  this  project  was  completed  at  the  beginning  of  this  last  contract 
year.  One  technical  paper  on  this  research  is  now  being  prepared  for 
publication. 

V.  AUTOMATION  OF  ARC  WELDING 

Automation  is  one  area  of  welding  technology  that  promises  to  pay 
tremendous  dividends  if  high  sensitivity,  rugged  sensors  can  be  developed. 

One  possible  sensor  is  the  arc  itself,  which  is  known  to  respond  to 
geometric  and  chemical  discontinuities  at  frequencies  up  to  20  kH^.  In 
one  study,  digital  signal  analysis  is  being  used  to  measure  these  responses. 

In  a  second  effort,  the  effect  of  convection  on  weld  pool  shape  and  defect 
formation  is  being  studied.  This  latter  topic  is  of  great  importance  in 
developing  processes  capable  of  producing  weld  pools  of  desired  size  and 
quality. 

A.  Signal  Analysis  as  a  Weld  Process  Monitor 

One  of  the  difficulties  associated  with  automation  of  arc  welding 
processes  is  the  lack  of  a  reliable  in-process  weld  quality  sensor.  The 
intense  heat,  light,  and  electric  fields  in  the  vicinity  of  the  arc  form  a 
hostile  environment  to  nearly  all  sensors.  By  using  the  arc  itself  as  a 
sensor,  this  difficulty  is  overcome.  The  purpose  of  this  research  project 
is  to  develop  a  weld  quality  sensor  based  on  the  electrical  properties  of 
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the  welding  arc.  Digital  signal  processing  of  the  arc  voltage  and  current 
is  an  effective  means  of  determining  these  properties.  During  previous 
contract  periods,  the  hardware  and  software  necessary  to  acquire  and  analyze 
the  signals  has  been  developed.  Preliminary  studies  have  shown  the  optimum 
current  waveform.  During  the  present  contract  period,  actual  measurements 
of  weld  properties  have  been  obtained. 

In  order  to  study  the  electrical  properties  of  the  arc,  experiments 
are  performed  with  welding  current  regulated  to  within  1%  of  the  selected 
value  and  a  constant  arc  length.  The  current  is  considered  to  be  the  input 
to  the  arc/puddle  system  and  the  arc  voltage  is  considered  to  be  the  output. 
These  voltage  and  current  signals  are  analyzed  to  determine  the  frequency 
response  function  of  the  arc  puddle  system  and  the  electrical  properties  of 
the  arc  are  determined  by  calculating  the  current-voltage  transfer  function 
of  the  system  (2) . 

The  experimental  work  performed  during  the  past  year  can  be  classified 
into  two  groups.  The  first  group  consists  of  experiments  designed  to  show 
the  effects  of  weld  defects  on  the  electrical  properties  of  the  arc,  and  the 
second  consists  of  experiments  designed  to  study  the  dynamics  of  the  weld 
pool.  The  current  waveform  and  sampling  parameters  have  been  optimized 
during  the  experimental  program  to  improve  the  signal-noise  ratio  of  the 
system  gain  function.  Thus,  the  quality  of  the  received  signals  is  now  very 
high,  providing  a  high  degree  of  confidence  in  the  results.  Although  the 
sensitivity  of  the  technique  remains  to  be  assessed,  the  capability  to 
measure  puddle  oscillations  has  been  demonstrated. 

A  limited  ability  of  this  analysis  technique  to  detect  changes  in  weld 
chemistry,  weld  geometry,  and  weld  size  has  been  shown.  Work  is  continuing 
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in  an  effort  to  identify  the  limits  of  sensitivity  and  practicality  of 
this  technique  as  a  means  of  sensing  weld  quality.  A  thesis  on  this 
topic  and  one  or  two  papers  will  be  completed  in  the  next  contract  year. 

B.  Convection  in  Arc  Weld  Pools 

Recent  work  by  many  investigators  is  emphasizing  the  importance  of 
fluid  convection  on  arc  weld  pool  geometry.  Related  work  supported  by 
the  U.S.  D.O.G.  in  this  laboratory  resulted  in  a  technical  paper  which 
recently  won  an  AWS  award  (3) . 

Work  on  this  project  has  concentrated  mainly  in  the  high  current 
range  (250-500A).  At  these  currents,  a  deep,  surface  depression  is  formed 
which  can  be  explained  by  a  rapid  vortex  flow  in  the  pool.  This  theory 
has  been  explained  in  a  previously  published  paper  (4) , 

One  technical  paper  is  scheduled  to  be  published  soon  on  arc  force 
measurements  (Appendix  H)  in  support  of  this  work.  This  paper  presents 
arc  force  measurements  up  to  600A  as  an  extension  of  a  previously 
published  technical  note  (5).  A  Thesis  on  the  entire  project  will  be 
completed  in  the  next  contract  year,  and  one  or  two  additional  papers 
are  being  prepared.  This  project  will  be  continued. 


SUMMARY 


During  the  past  year,  significant  progress  has  been  made  in  many 
areas  of  this  project.  The  research  conducted  over  this  time  involved 
spectroscopy  of  arc  welding,  welding  of  heavy  section  titanium,  fracture 
toughness  of  HY-80  weldments,  and  automation  of  arc  welding. 
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Metal  Vaporization  from  Weld  Pools 

A.  BLOCK-BOLTEN  and  T.  W.  EAGAR 

Experimental  studies  of  alloy  vaporization  from  aluminum  and  stainless  steel  weld  pools  have  been 
made  in  order  to  test  a  vaporization  model  based  on  thermodynamic  data  and  the  kinetic  theory  of 
gases.  It  is  shown  that  the  model  can  correctly  predict  the  dominant  metal  vapors  that  form  but  that 
the  absolute  rate  of  vaporization  is  not  known  due  to  insufficient  knowledge  of  the  surface  temperature 
distribution  and  subsequent  condensation  of  the  vapor  in  the  cooler  regions  of  the  metal.  Values  of  the 
net  evaporation  rates  for  different  alloys  have  been  measured  and  are  found  to  vary  by  two  orders  of 
magnitude.  Estimated  maximum  weld  pool  temperatures  based  upon  the  model  are  in  good  agreement 
with  previous  experimental  measurements  of  electron  beam  welds. 


I.  INTRODUCTION 

Loss  of  alloying  elements  from  die  weld  pool  due  to 
vaporization  is  important  for  a  number  of  reasons.  Firstly,  if 
the  loss  is  great  enough,  the  mechanical  properties  of  the 
weld  may  be  impaired.1  Secondly,  the  composition  of  a 
welding  arc  plasma  influences  the  temperature  of  the  arc,2 
arc  stability,  and  fume  formation.1  Thirdly,  it  has  been 
shown  that  vaporization  places  an  upper  limit  on  the  tem¬ 
perature  produced  on  the  surface  of  the  metal  due  to  evapo¬ 
rative  cooling.4,5 

A  previous  paper  has  presented  a  formalism  for  calcu¬ 
lation  of  partial  pressures  of  metal  vapors  above  steel  weld 
pools.5  This  analysis  provided  an  estimate  of  the  power  lost 
by  evaporation  as  well  as  an  upper  bound  on  the  surface 
temperature  of  steel  weld  pools  as  a  function  of  alloy  com¬ 
position.  In  the  present  paper,  this  analysis  is  extended  to 
evaporation  from  aluminum  and  copper  alloy  weld  pools 
where  different  metal  vapors  dominate.  The  results  of  the 
calculations  are  then  compared  with  experimental  results 
from  both  aluminum  and  steel  weld  metals. 


II.  EXPERIMENTAL  PROCEDURE 

Both  steel  samples  and  aluminum  alloy  samples  were 
subjected  to  extended  time  welding  in  a  specially  adapted 
chamber  which  was  coupled  to  a  direct  reading  emission 
spectrometer.  The  chamber  was  supplied  with  a  lens  guiding 
the  light  to  the  spectrometer  grating.  A  mirror  imaging  de¬ 
vice  was  provided  to  ensure  uniform  position  of  the  work 
and  the  electrode.  The  hearth  was  water  cooled  and  the 
entire  system  was  purged  with  argon  flow.  The  rotating 
water-cooled  copper  hearth  shown  in  Figure  1  was  em¬ 
ployed  in  the  case  of  12.S  cm  diameter  304  steel  samples, 
but  the  rotation  has  not  proven  useful.  The  slight  deviations 
of  the  rotating  sample  cause  periodic  sinusoidal  changes  in 
the  arc  length  and  corresponding  oscillations  in  the  spectral 
signal.  Therefore,  this  experiment  could  be  characterized 
only  by  lower  and  upper  limits  of  spectrographic  signal, 
voltage,  current,  and  final  weld  chemical  composition,  and 
not  by  a  uniquely  defined  value  of  each  of  those  parameters. 
As  a  result,  the  experimental  data  presented  in  this  paper 

A.  BLOCK-BOLTEN,  Research  Scientist,  is  now  with  the  Depart¬ 
ment  of  Metallurgy  and  Materials  Science,  University  of  Toronto. 
T.  W.  EAGAR.  Associate  Professor,  Materials  Engineering,  is  with  the 
Department  of  Materials  Science  and  Engineering,  Massachusetts  Insti¬ 
tute  of  Technology,  Cambridge,  MA  02139. 

Manuscript  submitted  September  23,  1983. 
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Fig.  I  — Water-cooled  turntable  and  gas-cooled  electrode  holder. 


were  obtained  from  stationary  arc  welds  where  the  arc 
length  was  more  controllable. 

The  0.35  to  1.10  g  steel  samples  consisted  of  1.6  mm 
diameter  wires  (502, 505, 5151,  and  5212  steels)  or  1. 1  mm 
wire  (410  steel)  or  2.4  mm  diameter  wires  (308L  and  309L 
steels).  The  thoriated  tungsten  electrodes  of  1.6  mm  di¬ 
ameter  were  mounted  as  shown  in  Figure  1.  Each  steel 
sample  was  weighed  before  and  after  each  experiment,  yet 
only  the  308L  and  309L  stainless  steel  samples  gave  reason¬ 
able  weight  loss  results  (Figure  2),  as  other  steels  were  too 
readily  oxidized  in  spite  of  the  welding  grade  argon  environ¬ 
ment.  The  small  amount  of  impurities  in  this  gas  stream  and 
the  extended  time  of  these  tests  created  a  significant  weight 
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Weight  loss  % 


Fig.  2 — Weight- loss  experiment.  I  — weight  Iom  curve  for  308L  stainless 
steel,  2 — for  309L  stainless  steel;  average  voltage  II  V,  current  set  at 
15  A;  avenge  sample  size  1  g. 


loss  or  gain  by  oxidation  over  the  course  of  the  experiment 
in  these  other  steels.  The  compositions  of  the  steel  and 
Al-alloy  samples  used  are  given  in  Table  I.  The  sizes  of  the 
Al-alloy  samples  are  specified  in  the  caption  of  Figure  8. 

Atomic  absorption  chemical  analysis  was  performed  on 
buttons  from  the  melted  steel  wires  and  on  shavings  drilled 
out  of  craters  of  the  Al-alloy  samples  after  melting. 

The  spectrographic  signals  were  ratioed  relative  to  the 
most  abundant  element  of  the  alloy:  Fe  in  steels  and  A1  in 
Al-alloys,  through  a  system  of  analog  dividers.  The  spectro¬ 
graphic  signals  from  die  photomultipliers  were  measured  by 
Keithley  Model  480  picoamperometers. 


III.  PREDICTION  OF 
VAPORIZATION  TENDENCY 

As  shown  in  a  previous  paper,3  the  logarithm  of  the  partial 
pressure  of  an  alloy  component  in  the  gas  phase  is  propor¬ 
tional  to  die  sum  of  the  logarithms  of  the  standard  pressure 
of  the  pure  element  and  the  activity  of  the  element  in  the 
alloy,  i.e. , 

log  p„  =  log  pi  +  log  aA  [1} 

where 

pA  is  the  partial  pressure  of  element  A  in  the  gas  phase. 
Pa  is  the  standard  pressure  of  pure  A,  and 
aA  is  the  activity  of  A  in  the  alloy. 

The  vaporization  rate  can  be  predicted  from  the  kinetic 
theory  of  gases,6 

rA  =  44.33  lpA^^  [g  •  s'1  •  cm’2]  12] 

where 

r*  is  the  rate  of  evaporation  of  element  A, 

Mk  is  the  molecular  weight  of  A, 


T  is  the  absolute  temperature,  and 
PA  has  units  of  atmospheres. 

The  evaporation  power  loss  PL  is  then 

Pi  =  rK(LK  -  AWa)  [watt  •  cm-2]  [3J 


where 

La  is  the  heat  of  evaporation  of  pure  A  and 

3/7a  is  the  partial  molar  heat  of  mixing  of  A  in  the  alloy. 

The  evaporative  energy  loss,  EL,  is 

El  =  rA(LA  ”  SffA)r[watt  cm'2  s],  [4] 

where  t  is  time.  In  most  cases  AWa  is  small  compared  to  the 
heat  of  evaporation.  If  the  element  of  interest  is  the  solvent 
such  as  iron  in  steel  or  aluminum  in  aluminum  alloys, 
can  be  neglected. 

Using  Eq.  [1],  it  is  possible  to  construct  pressure- 
temperature  diagrams  for  aluminum,  copper,  and  iron  base 
alloys  as  shown  in  Figures  3,  4,  and  3.  These  diagrams 


Fig.  3  —  Aluminum  alloys:  vapor  pressures  for  some  consliluent  alloy 
additions. 


Table  I.  Compositions  of  Starting  Materials  Used  in  This  Study 


Steels:  Compositions  As  Delivered 

At  Alloys:  Composition  As  Delivered 

Steel  Number 

Pet  Fe 

Pet  Mn 

PctCr 

Alloy  Number 

Pet  Zn 

Pet  Mg 

410 

86.70 

0.45 

12.03 

7075 

5.82 

2.30 

502 

94.87 

0.45 

4.00 

5083 

0.014 

4.06 

505 

89.70 

0.45 

9.06 

5456 

0.019 

5.20 

5151 

98.05 

0.54 

0.99/1.37 

6061 

0.063 

0.95 

5212 

97.30 

0.55 

1.91 

2024 

0.058 

1.46 

308L 

67.25 

1.91 

21.00 

1100 

<0.010 

<0.010 

309L 

59.65 

1.79 

26.08 

304 

68.80 

1.75 

18.66 
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Fig.  4 — Copper  alloys:  vapor  pressures  for  some  constituent  alloy 
additions. 


are  based  on  tabulated  values  of  activity  coefficients.12  To 
extrapolate  Hultgren’s  activity  (op)  data  to  other  ( a„ )  tem¬ 
peratures,  log  a„  =  log  Oo  +  AG"  (To  -  7n)/4.575  ToTn, 
AG"  was  assumed  identical  to  A//  and  AG"  or  AT/  was 
used,  whichever  was  available.  AG"  is  partial  excess  free 
energy  of  mixing. 

As  shown  previously,  mild  interaction  effects  between 
the  alloying  elements  in  ternary  and  higher  order  systems 
can  often  be  neglected;5  however,  strong  interactions  where 
compounds  tend  to  form  in  the  liquid  metal  cannot  be 
neglected. 

As  shown  previously  for  steels,5  Eq.  [4]  can  be  com¬ 
bined  with  the  arc  surface  energy  distributions  of  Nestor7 
to  produce  upper  limits  on  the  surface  temperature  of  the 
weld  pool.  Figure  6  shows  such  an  upper  bound  analysis  for 
a  number  of  aluminum  alloys.  The  upper  bound  analysis 
assumes  that  all  of  the  power  input  to  the  weld  pool  is  lost 
by  vaporization  of  the  metal.  Since  this  is  clearly  not  true, 
a  least  upper  bound  on  the  temperature  of  the  weld  pool 
surface  can  be  obtained  if  one  knows  what  fraction  of  the 
total  input  power  is  lost  by  evaporation.  As  will  be  shown 
later,  it  is  difficult  to  measure  this  quantity  with  precision. 
The  concept  of  the  upper  bound  temperature  of  the  pool  is 
outlined  in  a  previous  paper.5 


IV.  CORRELATION  OF  THEORY 
WITH  EXPERIMENT 

As  noted  in  the  experimental  procedure  section,  a  number 
of  steel  and  aluminum  alloys  were  arc  melted  for  unusually 
long  times  in  order  to  evaporate  enough  alloying  compo¬ 
nents  to  be  detectable  by  chemical  analysis.  Due  to  prob¬ 
lems  of  oxidation,  only  two  of  the  stainless  steel  samples 
gave  consistent  results.  The  Mn  and  Fe  spectrograph ic 
signals  normalized  to  the  relatively  stable  Cr  signals  as  a 


13 


Fig.  5  —  Steels:  vapor  pressures  for  some  constituent  alloy  additions. 
Numbers  on  curves  represent  weight  percentages. 
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function  of  arc  melting  time  are  shown  in  Figure  7.  It  will 
be  noted  that  Mn  is  the  primary  element  lost  as  predicted 
by  Figure  5. 

In  order  to  test  the  model  more  completely,  a  series  of  six 
aluminum  alloys  with  varying  zinc  and  magnesium  contents 
were  selected.  Aluminum  alloys  were  selected  because  the 
Zn  and  Mg  have  very  high  vapor  pressures  as  compared  to 
aluminum;  and  as  seen  in  Figure  6,  this  significantly  lowers 
the  upper  bound  temperature. 

The  aluminum  alloys  can  be  separated  into  four  groups, 
where 

I  Zn  vapor  dominates, 

II  Mg  vapor  dominates, 

III  nearly  equivalent  Zn  and  Mg  vapor  pressure  exists, 

IV  A1  vapor  dominates. 

Alloy  7075  belongs  in  group  1  as  seen  in  Figure  8.  If 
we  compare  Table  II  and  Figure  3,  we  see  that  the  high 
zinc  concentration  in  7075  alloy  places  the  zinc  isopleth  of 
Figure  2  well  above  the  Mg  isopleth  for  this  alloy. 

Alloys  5083  and  5456  belong  to  group  II  as  seen  in 
Figures  9  and  10.  This  is  again  consistent  with  the  predic¬ 
tions  made  from  Table  II  and  Figure  3. 

Alloys  6061  and  2024  belong  to  group  III  as  shown  in 
Figures  11  and  12,  respectively.  This  result  is  not  as  easily 
predicted  from  Figure  3  since  the  exact  composition  of  the 
alloy  and  the  temperature  of  the  weld  pool  surface  are  not 
known  precisely.  Nonetheless,  the  spectrographic  data  of 
Figure  13  clearly  show  a  simultaneous  rise  in  the  spectral 
signals  for  both  Mg  and  Zn.  It  is  interesting  to  note  that 
Zn  typically  reaches  its  maximum  presence  in  the  vapor  later 
than  Mg,  probably  due  to  kinetic  limitations  since  Zn  is 
more  dilute  in  the  alloy.  The  gradual  rise  of  both  elements 


Fig.  7  —  Decay  of  normalized  signal  ratios  of  iron  to  chromium  and  manganese  to  chromium 
(Cr  content  does  not  change  significantly)  for  three  different  steels  As  can  be  seen,  Mn  domi¬ 
nates  over  the  iron  in  the  vapor  phase  Only  after  long  periods  of  time  does  iron  start  to 
dominate  in  the  SI5I  and  $02  steels. 


NM- TEMPERATURE  DISTRIBUTION  CSM-TEMPERATURE  DISTRIBUTION 


Fig.  6 — Normal  mode  and  cathode  spot  mode  temperature  distributions  for 
the  5456  Al-alloy.  The  temperature  is  plotted  vs  non-dimensional  distance 
r/p  from  the  center  of  the  arc;  p  is  defined  as  in  Ref.  5  as  p  =  VV/2 irC 
where  V  is  the  total  power  of  the  arc  (in  watts)  and  C  is  the  peak  power 
density  in  watt  •  cm"2.  The  plots  are  based  on  Fig.  3  and  on  Nestor’s7 
power  distribution  curves.  Should  100  pet  of  energy  go  into  evaporation  of 
pure  aluminum,  point  A  on  temperature  ordinate  would  correspond  to  the 
expected  maximum  temperature  of  the  weld  pool.  Points  B  and  C  mark 
10  pet  and  1  pet  energy  going  into  evaporation,  respectively.  Curves  I  are 
normal  distribution  (of  energy)  curves  for  pure  aluminum.  Curves  2  are 
those  based  on  Nestor  for  pure  Al.  Curves  3  are  for  evaporation  of  0.25  pet 
Zn  in  Al,  and  curves  4  and  5  are  for  evaporation  of  4.7  pet  Mg  in  Al  and 
5.5  pet  Mg  in  Al,  respectively. 
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Fig.  8 — Loss  of  Zn  vs  energy  input  in  [KVA  min]  from  aluminum  alloy  7075.  No  change  in  Mg  was  found. 
Original  contents  of  Mg  and  Zn  were  2.3  pet  and  5.6  pet,  respectively.  The  voltage  during  tests  varied  be¬ 
tween  12.5  and  15  V;  current  was  set  between  125  A  and  ISO  A;  the  longest  experiment  lasted  45  min.  The 
size  of  die  samples  was  52  x  51  x  12.7  mm.  Rates  of  evaporation  for  this  and  other  alloys  are  given 
in  Table  III. 

ible  II.  Calculated  Center  Line  Surface  Temperature  of  WeM  Pools  for  Six  Aluminum  Alloys 


1 100  Alloy.  0  to  0. 1  *  Pet  Zn,  Bal.  A1 


Power  Vapor  t  °C  Range 

Density,  Pet  Domination  NM/CSM 


100  A1  1935*/2091 

10  Zn*.  A1  1364*/1614 

1  Zn*,  A1  968*/1530 

5456  Alloy,  0.25  Pet  Zn; 

4.7  to  5.5*  Pet  Mg 


Power  Vapor  t  °C  Range 

Density,  Pet  Domination  NM/CSM 


100  Mg  1220*/ 14-51  100  Mg  1513 

10  Mg  964*/1104  10  Zn,  Mg*  1154 

1  Mg  783*/  891  1  Zn  870 


2024  Alloy,  0.25  Pet  Zn; 
1.2  to  1.8*  Pet  Mg 


Power 
Density.  Pet 


100 

10 

1 


Vapor 

Domination 


Mg 

Mg 

Zn,  Mg*CSM 


t  °C  Range 
NM/CSM 


1422*/ 1751 
1101*/1290 
870*/  975 


6061  Alloy,  0.25  Pet  Zn; 
0.8  to  1.2*  Pet  Mg 


Power 
Density,  Pet 


100 

10 

1 


Vapor 

Domination 


Mg 

Zn,  Mg* 
Zn 


Normal  Mode  NM  200  A,  14.0  V,  C  =  2180  [wan  •  cm'!l 

Cathode  Spot  Mode  CSM  200  A,  14.8  V,  C  =  6120  [wan  ■  cm'J]  peak  power  density 
•Note  that  the  lower  calculated  temperatures  correspond  to  the  higher  alloy  compositions. 


I  °C  Range 
NM/CSM 


1513*/1882 
1154*/ 1358 
870*/  994 


5083  Alloy,  0.25  Pet  Zn; 
4  to  4.9*  Pet  Mg 


Power  Vapor  t  °C  Range 

Density,  Pet  Domination  NM/CSM 


100  Mg  1268*/1491 

10  Mg  994*/ll26 

1  Mg  807*/  905 

7075  Alloy,  5.1  to  6.1*  Pet  Zn; 

2.1  to  2.9  Pet  Mg 


Power  Vapor  t  °C  Range 

Density,  Pet  Domination  NM/CSM 


100  Zn  1057*/ 1242 

10  Zn  795*/  909 

1  Zn  612*/  695 


! 


VIt[KVAmin] 


Fig.  10 — Aluminum  alloy  5456.  Loss  of  Mg  vs  energy  input.  The  changes 
in  Zn  cootent  occur  only  in  the  first  seconds  of  the  experiment,  then  Zn 
remains  constant.  Original  Mg  content  5.8  pet.  Avenge  voltage  was 
15  V,  and  current  was  set  between  125  and  150  A.  The  longest  experiment 
lasted  45  min. 


is  most  likely  due  to  increase  in  the  size  of  the  weld  pool 
with  time.  Alloy  1100  is  the  only  sample  tested  which  be¬ 
longs  to  group  IV  of  the  series. 

Summarizing,  Figures  8  to  12  show  the  analyses  taken 
from  die  weld  pools  of  different  Al-alloys  during  extended 
time  welding.  The  707S  alloy  with  Zn-vapor  domination 
(Figure  8)  exhibits  changes  only  in  Zn  content,  while  Mg 
content  stays  constant.  In  Figures  9  and  10,  the  Mg-vapor 
dominated  50S3  and  S456  alloys  show  Mg-content  falling 
off,  while  Zn -con  tent  stays  constant  or  below  the  detection 
level.  Figures  11  and  12  for  alloys  6061  and  2024  prove 
that  both  these  alloys  exhibit  mixed  vapor  domination  and 
both  Zn  and  Mg  contents  change  with  time.  In  addition, 
SEM  semiquantitative  analyses  have  also  shown  changes 
in  Zn  only  for  the  707S  alloy,  changes  in  Mg  only  for  5083 
and  5456  alloy,  and  changes  both  in  Mg  and  Zn  content  for 
the  2024  and  6061  alloys. 


Fig.  II — Aluminum  alloy  6061.  Lou  of  Mg  and  Zn  vi  energy  input.  Original  Mg  content  was  0.95  pet 
and  the  original  Zn  content  0.06  pet.  Average  voltage  14  V;  the  current  was  set  at  150  A.  Longest  experiment 
lasted  45  min. 


Table  III.  Net  Energy  Lost  from  Weld  Pool  by  Evaporation  of  Alloying  Elements 


Group 

Al-Alloy 

Evaporating 

Element 

Apparent 
Initial  Net 
Evaporation 
Rate 

(Mgs-1  cm'2] 

Apparent  Long- 
Term  (Mean) 
Evaporation 
Rate 

[#xgs~'  cm"2] 

Theoretical 

Evaporation 

Energy 

(t-Atf)" 

U*r'i 

Apparent 
Evaporation 
per  Gross 
Energy  Input* 
iMg  J”  x  H)3) 

Percent  of 
Arc  Energy 
Lost  by 
Evaporation 

Source 
of  Data 
(Figure  No.) 

I 

Zn 

94.9 

30.5 

631.1 

91.4 

0.0144 

8 

Mg 

— 

— 

173.5 

— 

— 

II 

Zn 

(11.3) 

631.1 

9 

Mg 

5.6 

12.1 

173.5 

5456 

Zn 

— 

631.1 

— 

— 

Mg 

9.7 

12.7 

173.5 

49.0 

Al** 

92.8 

— 

— 

III 

Zn 

1.2 

1.2 

631.1 

3.6 

11 

Mg 

1.2 

173.5 

3.6 

Zn 

mBBm 

1.3 

631.1 

3.8 

12 

Mg 

Wmmt 

3.4 

173.5 

3.4 

*grou  energy  input.  Refs.  8.  9 
••from  soot  formation  and  Al-ball  growth 
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F>g  12 —  Aluminum  alloy  2024.  Loss  of  Mg  and  Zn  content  vs  energy  input.  Original  Mg  content  1.5  pet, 
Zn  content  0.06  pet.  Average  voltage  was  14  V,  and  the  current  was  set  at  ISO  A.  Longest  experiment 
lasted  30.2  min. 


6061  Al  Alloy 


Fig.  13  —  Signal  ratios  of  Mg  and  Zn  normalized  to  aluminum,  as  typically 
observed  for  the  606 1  aluminum  alloy  The  Mg/AI  signal  ratio  is  magnified 
10  times. 


From  Figures  8,  9,  10,  11,  and  12,  values  of  apparent 
initial  evaporation  rate,  apparent  long  term  (mean)  evapo¬ 
ration  rate,  and  apparent  evaporation  per  unit  of  energy 
have  been  calculated.  The  latter  value  is  normalized  to  gross 
energy  input'9  which  is  half  of  the  total  arc  power  input, 
as  the  other  half  is  lost  to  the  environment  and  does  not 
enter  the  weld  pool.  All  these  values,  together  with  percen¬ 
tages  (referred  to  theoretical  evaporation  per  unit  energy 
(L  -  57/r'[r'l)  are  displayed  in  Table  III.  It  is  seen  that 
the  net  power  lost  by  vaporization  is  much  less  than  one 
percent  of  the  total  power  input  to  the  metal.  As  will  be 
discussed  later,  this  is  much  less  than  the  absolute  evapora¬ 
tive  power  loss;  however,  much  of  the  absolute  evaporative 
power  is  regained  by  condensation  of  the  vapors  in  cooler 
regions  of  the  weld  pool. 
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V.  DISCUSSION 

The  preceding  results  indicate  that  the  model  can  effec¬ 
tively  use  existing  thermodynamic  data  to  predict  the  domi¬ 
nant  metal  vapors  above  the  weld  pool;  however,  it  would 
be  useful  if  one  could  predict  weld  pool  surface  temperature 
and  hence  the  rate  of  alloy  element  loss.  This  proves  to  be 
a  much  more  difficult  task. 

Quigley,  el  al. 9  have  estimated  the  evaporative  power 
loss  from  a  gas  tungsten  arc  to  be  of  the  order  of  2  pet  of  the 
total  power.  If  this  is  true,  it  would  reduce  the  upper  bound 
temperatures  based  on  100  pet  evaporative  power  loss  as 
shown  in  Figure  6  and  in  Table  II.  Indeed,  if  one  were  to 
assume  absolute  evaporative  power  loss  to  lie  between  1  and 
10  pet  of  the  total  power,  the  surface  temperature  limit 
would  drop  by  200  to  500  °C  for  aluminum  alloys.  Table  II 
shows  that  the  surface  temperature  of  aluminum  alloy  arc 
welds  could  vary  from  just  above  the  melting  temperature  to 
700  °C  superheat  depending  on  alloy  composition  and  the 
amount  of  heat  lost  by  vaporization. 

As  shown  in  Table  III  and  Figures  8  though  12,  weight 
loss  measurements  based  on  alloying  element  loss  from  the 
weld  samples  gives  net  evaporative  losses  between  zero  and 
95[jtg  *  s'1  •  cm-2].  This  corresponds  to  maximum  net 
vaporization  losses  of  roughly  0.015  pet  of  the  net  incident 
power;  however,  it  is  believed  that  more  metal  than  this  is 
vaporized  from  the  center  of  the  pool  and  then  recondenses 
on  cooler  regions,  thus  recycling  through  the  system.  This 
is  indirectly  confirmed  by  spectrographic  monitoring  of  cal¬ 
cium  vapor  in  304  stainless  steels.  The  calcium  spectro¬ 
graphic  lines  remain  strong  and  do  not  diminish  after 
15  minutes  of  welding  even  though  there  is  only  10  ppm  Ca 
in  the  steel.  The  only  way  this  small  amount  of  Ca  could 
remain  after  such  a  long  time  is  if  it  is  “recycled”  into  the 
metal  by  condensation.  Although  it  is  not  possible  to  mea¬ 
sure  how  much  vapor  condenses  back  into  the  pool  rather 
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than  is  lost,  it  is  believed  that  a  larger  fraction  of  Zn  recon¬ 
denses  back  to  the  pool  than  is  lost  from  the  system.  In  a 
sense,  this  is  confirmed  by  the  temperatures  calculated  for 
the  7075  alloy  in  Table  II.  If  the  temperature  were  low 
enough  for  only  0.015  pet  of  the  power  to  be  lost  by  vapor¬ 
ization,  then  the  liquid  would  be  cool  enough  to  solidify. 
It  is  clear  that  at  least  several  percent  of  the  power  at  the 
center  of  the  7075  alloy  pool  must  be  lost  by  zinc  vapori¬ 
zation.  Much  of  this  power  is  recovered  at  the  edges  of  the 
pool  as  the  vapor  condenses.  Because  of  this  evaporation- 
condensation  process  it  is  very  difficult  to  determine  the 
absolute  vaporization  rate  with  any  accuracy. 

Schauer  et  al.  used  an  infrared  pyrometer  to  measure  the 
weld  pool  surface  temperature  of  electron  beam  welds  in 
steel  and  aluminum."  If  one  assumes  that  the  power  lost 
from  the  electron  beam  weld  pool  is  6000  watts  *  cm"2,  the 
measured  temperatures  are  in  reasonable  agreement  with  the 
values  predicted  by  the  model  presented  here  as  seen  in 
Table  IV.  Indeed,  these  results  are  perhaps  the  most  conclu¬ 
sive  in  showing  that  evaporative  power  losses  set  an  upper 
limit  on  the  surface  temperature  of  a  weld  pool.  In  general, 
these  power  losses  probably  lie  between  1  and  10  pet  of  the 
total  power,  but  the  value  varies  somewhat  with  alloy  con¬ 
tent  and  with  input  power  density.  These  losses  provide  a 
prediction  of  maximum  arc  weld  pool  surface  temperatures 
between  the  melting  point  and  1600  °C  for  aluminum  alloys 
and  between  2000  °C  and  2500  °C  for  steels.  The  high  en¬ 
ergy  density  processes  such  as  laser  and  electron  beam  can 
produce  even  higher  surface  temperatures  with  greater  rates 
of  alloy  vaporization  than  arc  welds. 


A  thermodynamic  model  of  vaporization  from  weld  pools 
has  been  presented  which  can  predict  the  dominant  metal 
vapors  that  form  during  welding  in  the  absence  of  strong 
compound  forming  elements  such  as  oxygen  or  nitrogen. 
The  evaporative  power  loss  in  arc  welding  lies  between 
1  and  10  pet  of  the  incident  power,  although  more  precise 
measurement  is  difficult  due  to  condensation  of  metal  va¬ 
pors  in  cooler  regions  of  the  weld  pool.  These  vaporization 
rates  place  upper  limits  on  the  surface  temperature  of  the 
weld  pool.  These  limits  are  in  good  agreement  with  mea¬ 
surements  on  different  alloys  when  using  electron  beam 
welding,  but  the  results  for  arc  welding  are  somewhat  lower 
than  the  temperatures  often  assumed  for  arc  weld  pools.  It 
is  concluded  that  the  metal  vaporization  from  the  weld  pool 
places  a  limit  on  the  maximum  temperature  of  the  pool 
which  is  significantly  less  than  the  boiling  temperature  of 
the  metal.  The  presence  of  volatile  alloying  elements  may 
further  reduce  this  maximum  temperature  limit. 

Table  IV  indicates  that  apparent  evaporation  rates  (initial 
and  mean)  confirm  Zn  domination  of  the  group  I  alloys.  Mg 
domination  in  group  II,  and  co-domination  of  Zn  and  Mg  in 
group  III  alloys. 

The  largest  absolute  evaporation  seen  occurs  in  the  Zn 
dominated  group  1  followed  by  Mg  dominated  group  II.  The 
apparent  percentages  of  the  total  energy  are  very  small  indi¬ 
cating  that  the  vapors  are  subjected  to  multiple  recycling 
during  extended  time  arc  welding. 


Table  IV.  Comparison  o t  Predicted  Welding  Temperatures  and  Temperatures  Measured  by  Schauer  et  al. " 


Al-Alloy,  Steel 
or  Metal 

Maximum  Predicted  Temperatures  (CSM) 
if  10  Pet  or  100  Pet  Energy  Went 
into  Evaporation 

Dominating  Vapor 
Species  at  CSM  if 

Schauer 's"  Electron 
Beam  Temperatures 

t  °C 

10  Pet  100  Pet 

Went  into  Evaporation 

Aluminum 

1775 

A) 

2091 

Al 

1100 

1614 

Zn 

1900  ±  100 

2091 

Al 

(1990  to  2055)* 

2024 

1290 

Mg 

1700  ±  100 

1751 

Mg 

5083 

1126 

Mg 

1250  ±  100 

1491 

Mg 

5456 

1104 

Mg 

1451 

Mg 

6061 

1358 

Mg 

1882 

Mg 

(1800  to  1890)* 

7075 

909 

Zn 

1242 

Zn 

(1380  to  1485)* 

Steels  with  0.5  pet  Mn 

2190 

Mn 

2290  ±  60;  HY-130st 

2520 

Fe 

Steels  with  2.0  pet  Mn 

2010 

Mn 

2100  ±  50;  304  st 

2520 

Fe 

20-6-9  steel 

Mn 

1820  ±  40 

Mn 

Tantalum 

Ta 

4400  t  150 

Ta 

•Data  from  Schaucr's  figures,  not  matching  tabulated  temperatures 
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ABSTRACT 


Metal  vapors  in  gas  tnngsten  welding  arcs  were  studied  in  order 
to  determine  the  effects  of  these  vapors  on  arc  properties  and  subse¬ 
quently  on  weld  bead  configuration.  Emission  spectroscopy  and 
monochromatic  photography  were  used  to  determine  the  composition  and 
distribution  of  vapors  in  arcs  on  stainless  steel.  It  was  found  that, 
in  addition  to  the  expected  species  (Fe,  Mn,  Cr) ,  calcium  and  aluminum 
vapors  were  also  detectable,  even  at  very  low  base-metal  concentra¬ 
tions.  All  metal  vapors  were  determined  to  be  concentrated  just  above 
the  weld  pool  in  high  current  arcs,  but  also  present  in  the  upper 
regions  of  low  current  arcs.  Aluminum  and  calcium  were  also  found  to 
vaporize  from  the  tungsten  electrode. 

Theoretical  calculations  of  the  plasma  transport  properties 
revealed  that  small  amounts  of  low  ionization  potential  elements  such 
as  aluminum  or  calcium  do  not  have  as  great  an  effect  on  the  electrical 
and  thermal  conductivities  as  has  been  previously  reported  if  the  pres¬ 
ence  of  other  metal  vapors  such  as  iron  or  manganese  is  also 
considered. 

It  is  therefore  concluded  that  the  effects  of  minor  elements  on 
arc  properties  nay  be  less  important  than  has  previously  been  believed 
in  explaining  the  variable  penetration  often  associated  with  minor  ele¬ 
ment  additions  to  the  base  metal,  and  that  weld  pool  convection  effects 
such  as  surface  tension  modifications  are  probably  more  important. 
However,  the  effects  of  vapors  emitted  by  the  tungsten  electrode  may 
have  a  great  effect  on  arc  properties,  as  the  shielding  gas  is  other¬ 
wise  free  of  contaminants  in  the  upper  regions  of  the  arc. 


Thesis  supervisor:  Thomas  W.  Eagar 

Title:  Associate  Professor  of  Materials  Engineering 


Appendix  C 

(presented  at  AIME  Annual  Meeting,  New  York,  NY,  February  1985) 

CHARACTERIZATION  OF  MOLTEN  WELD  POOL 
by 

J.  Szekely,  T.  W.  Eagar,  J.  McKelliget,  and  G.  Dunn 

ABSTRACT 

A  complete  understanding  of  heat  and  fluid  flow  in  the  molten  weld 
pool  requires  an  understanding  of  processes  occurring  within  the  arc  plasma 
itself.  Through  the  simultaneous  solution  of  Maxwell's  equations,  the 
Navier-Stokes  equation  and  the  heat  balance  equation  in  the  plasma,  one 
can  predict  the  temperature,  fluid  flow  and  electric  current  fields  in  the 
plasma.  This  permits  one  to  estimate  the  fraction  of  heat  conducted  to  the 
workpiece  by  the  current  and  by  convection.  It  is  further  shown  that  these 
fields  can  be  perturbed  by  the  introduction  of  metal  vapors  within  the  plasma. 
Experimental  studies  confirm  that  these  vapors  are  not  homogeneously  distri¬ 
buted  throughout  the  plasma,  hence,  the  complete  solution  to  this  problem 
becomes  formidable  although  not  intractable. 
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SPECTROMETRIC  TEMPERATURE  MEASUREMENTS 
OF  THERMAL  RADIATORS 
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on  May  4,  1984  in  partial  fulfillment  of  tbe  requirements 
for  tbe  Degree  of  Doctor  of  Philosophy  in  Metallurgy 


ABSTRACT 


It  is  demonstrated  that  the  temperature  of  a  thermal  radiator  can  be 
determined  by  multiple  spectral  radiance  measurements  without  prior 
knowledge  of  the  emissiwity  of  the  source.  This  new  passive  measurement 
technique  assumes  only  that  a  smooth  function  exists  between  spectral 
emissivity  and  wavelength.  The  spectral  radiance  values  are  fit  to  a 
Planck  radiation  law  relation  to  yield  the  temperature  of  the  source. 
Computer  simulations  are  included  to  show  that  relative  errors  in  the 
tea^erature  measurements  are  generally  an  order  of  magnitude  less  than 
relative  errors  in  the  spectral  radiance  measurements  and  in  the  simul¬ 
taneously  calculated  spectral  emissivity  values. 

A  Multichannel  Infrared-Red  Temperature  Micro-Analyzer  (MIRTMA) 
system  is  described.  This  early  prototype  instrument  is  capable  of  moni¬ 
toring  temperatures  above  1100  K  with  a  spatial  resolution  of  100  pm  and 
a  temporal  resolution  of  3  s  using  spectral  radiance  measurements  at 
approximately  200  wavelengths  in  a  range  of  0.6-0. 8  pm.  Demonstrations 
of  this  equipment  on  a  heated  platinum  strip  source  are  described  and 
discussed.  The  temperature  measurements  of  these  sources  are  generally 
within  5%  of  the  actual  temperature,  but  can  be  within  1%  using  certain 
techniques . 

Capabilities  and  limitations  of  the  method  snd  the  MIRTMA  equipment 
are  presented.  The  various  calibrations  used  with  the  technique  are  also 
explained,  particularly  with  respect  to  corrections  for  the  nonlinear 
response  of  the  Silicon  Intensified  Target  detector  employed  as  part  of 
the  MIRTMA.  Potential  improvements  to  the  method  snd  instrument  based  on 
this  work  are  presented. 
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An  Improved  method  of  multi-wavelenth  pyrometry 

Cordon  B.  Hunter,  Cherly  D.  Allemand,  and  Thomas  W.  Eagar 

Materials  Processing  Center,  Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 


Abstract 

It  is  demonstrated  that  the  temperature  of  a  thermal  radiator  can  be  determined  by  curve 
fitting  techniques  using  multiple  spectral  radiance  measurements  without  prior  knowledge  of 
the  emissivity  of  the  source.  This  new  passive  measurement  technique  assumes  only  that  a 
smooth  function  exists  between  spectral  emissivity  and  wavelength.  The  spectral  radiance 
values  are  fitted  to  a  Planck  radiation  law  relation  to  yield  the  temperature  of  the  source. 
Error  analysis  shows  that  relative  errors  in  the  temperature  measurements  are  generally  an 
order  of  magnitude  less  than  in  the  spectral  radiance  measurements  and  in  the  simultaneous¬ 
ly  calculated  spectral  emissivity  values.  Computer  simulations  are  included  which  show  the 
effects  of  varying  different  parameters,  such  as  the  number  of  data  pairs,  the  wavelength 
range,  the  spectral  emissivity  behavior,  the  source  temperature,  and  the  measurement  noise, 
on  the  accuracy  of  the  temperature  determination.  Experimental  confirmation  of  this  tech¬ 
nique  is  presented  showing  temperature  measurements  within  1%  of  the  actual  temperatures  on 
a  platinum  surface  within  a  temperature  range  of  1273-1724  K. 

Introduction 


One  of  the  greater  problems  faced  by  spectrometric  methods  of  temperature  measurement  is 
characterization  of  the  emissivity  of  the  source.  The  most  widely  used  types  of  radiation 
pyrometry  are  based  on  an  approximation  of  the  Planck  radiation  law  and  an  assumption  of 
constant  spectral  emissivity  in  the  observed  wavelength  band.  Ratio  pyrometry  techniques, 
such  as  the  two  color  method,  are  examples  of  this1-  .  Other  methods  require  that  emissiv¬ 
ity  be  determined  separately  from  the  temperature  to  obtain  accurate  thermal  measurements. 

One  class  of  techniques,  referred  to  as  multi-wavelength  pyrometry,  has  been  developed  to 
calculate  both  the  temperature  and  emissivity  of  a  thermal  radiator  employing  one  set  of 
radiation  measurements  from  several  wavelengths3-6.  Typically,  these  techniques  involve 
assuming  a  functional  dependence  between  spectral  emissivity  and  wavelength,  and  then  using 
the  radiation  measurement  data  to  determine  the  adjustable  emissivity  parameters  and  tem¬ 
perature  by  solving  simultaneous  equations.  The  main  disadvantage  with  these  methods  is 
the  sensitivity  of  the  accuracy  to  radiation  measurement  errors*. 

An  improvement  to  this  approach  is  to  use  a  much  larger  number  of  radiation  measurements 
to  determine  the  temperature  and  emissivity  parameters  using  curve  fitting  techniques.  The 
sensitivity  to  measurement  errors  is  greatly  reduced  by  the  statistics  involved  with  curve 
fitting  a  large  number  of  measurements.  The  only  assumption  necessary  for  this  technique 
is  that  a  smooth  function  exists  between  spectral  emissivity  and  wavelength  over  the  wave¬ 
length  range  of  the  radiation  data. 


Theory 

Provided  that  the  source  is  a  thermal  radiator,  the  spectral  radiance  (N^  in  units  of 
W/cm2-um-sr)  emitted  is  function  of  the  wavelength  (A  in  urn)  of  the  radiation,  the  tempera¬ 
ture  (T  in  K)  of  the  source,  and  the  spectral  emissivity  of  the  emitting  surface  ac¬ 

cording  to  the  Planck  radiation  law: 

Nx  -  ex  C,  /(A5  (eC*/AT-l))  (1) 

where  C,  (11910  W-umVcm* -sr)  and  C2  (14387.9  ym-K)  are  the  Planck  radiation  constants1-2'7. 
Typical  Planck  radiation  law  curves  are  shown  ir.  Figure  1.  Emissivity  is  a  function  of 
many  variables  ,  however  all  but  wavelength  are  constant  if  the  spectral  radiance  measure¬ 
ments  are  simultaneous  and  optically  identical. 

Because  it  has  been  difficult  to  establish  an  analytical  expression,  spectral  emissivity 
has  been  represented  in  many  studies  by  a  simple  series  of  terms  with  adjustable  parameters, 
a  polynomial  for  instance,  valid  in  a  restricted  wavelength  range3'6.  Substituting  the 
assumed  spectral  emissivity  function  into  Equation  1  yields  the  spectral  radiance  relation 
as  a  function  of  only  wavelength  at  a  constant  temperature.  The  temperature  and  coeffi¬ 
cients  of  the  emissivity  relation  are  obtained  simultaneously  by  curve  fitting  the  spectral 
radiance  measurements  for  a  large  number  of  wavelengths  to  the  spectral  radiance  function. 
This  eliminates  the  need  to  determine  emissivity  separately  from  the  temperature. 
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Figure  1.  Planck  radiation  law  curves  with  a  constant  spectral  emissivity  (e-v  *  1)  for  four 
temperatures.  The  dotted  line  denotes  the  locus  of  curve  maxima  for  all  temperatures  and 
this  spectral  emissivity. 

The  insensitivity  of  the  temperature  determination  to  errors  in  the  spectral  radiance 
measurements  can  be  shown  analytically  by  differentiating  Equation  1  and  substituting  the 
relations: 

Z  =  C2  /  X  T 

Y  =  Z  /  (1  -  e~z) 
into  the  result  to  obtain: 

dN^/N^  *  dcA/cA  +  Y  dT/T  (4) 

For  the  range  of  wavelengths  and  temperatures  examined  in  this  study,  Y  is  about  one  order 
of  magnitude.  Thus,  a  relative  error  in  the  spectral  radiance  data  will  cause  at  most  an 
equal  relative  error  in  spectral  emissivity  or  about  an  order  of  magnitude  less  relative 
error  in  temperature.  The  precision  in  the  emissivity  and  temperature  determinations  are 
further  increased  due  to  the  statistical  nature  of  the  curve  fitting  analysis. 

As  can  be  deduced  from  Equation  4,  correlation  effects  between  the  adjustable  parameters, 
even  in  the  absence  of  experimental  error  in  the  spectral  radiance  (dN^/N^  *  0)  ,  can  cause 
errors  in  the  temperature  determination: 

deA/e^  -  -Y  dT/T  (5) 

Application 

The  eventual  goal  of  this  work  is  to  produce  an  instrument  and  technique  to  measure  sur¬ 
face  temperatures  of  molten  gas-tungsten  arc  weld  pools.  Most  previous  studies  of  weld 
pool  temperatures  have  employed  thermocoup les3“ 1  * .  These  techniques  have  been  far  from  op¬ 
timal  since  the  thermocouples  impede  fluid  motion  in  the  pool  and  may  alter  the  thermal 
profile'.  To  avoid  contact  with  the  weld  pool  some  investigators  have  used  total  radiation 
pyrometry  to  map  isotherms  on  and  around  the  weld  pool11”13,  however  the  difficulty  in 
determining  the  emissivity  has  prevented  them  from  accurately  associating  these  isotherms 
with  absolute  temperatures. 

Recent  advances  in  multichannel  photoelectric  detectors  allow  the  simultaneous  spectral 
radiance  measurements  at  multiple  wavelengths.  Coupled  with  the  rapid  data  processing 
capabilities  of  computers,  these  detectors  make  practical  the  improved  pyrometric  technique 
described  above.  This  technique  has  the  potential  for  eliminating  the  problems  encountered 
with  thermocouples  or  emissivity  in  other  studies,  as  well  as  making  accurate  thermal 
measurements. 

An  instrument  was  built  to  experimentally  demonstrate  the  technique  and  serve  as  a  proto¬ 
type  for  eventually  measuring  weld  pool  surface  temperatures14.  The  instrument  is  named  the 
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(2) 

(3) 


Multichannel  Infrared-Red  Temperature  Micro-Analyzer  (MIRTMA) .  The  MIRTMA  ia  divided  func¬ 
tionally  into  three  parts:  the  infrared  microscope,  the  detector  system,  and  the  data  pro¬ 
cessor.  The  infrared  microscope  is  comprised  of  the  spectrograph,  which  disperses  the 
radiation  onto  the  photodetector  surface,  and  of  the  focusing  optics,  which  image  a  100  wm 
effective  source  area  onto  the  spectograph.  The  detector  system  measures  the  photon  irra- 
diance  and  provides  experiment  sychronizati.on  capabilities,  while  the  data  processor  is  used 
to  calibrate  the  detector  response  to  spectral  radiance  and  analyze  the  calibrated  data  to 
determine  the  source  temperature.  The  detector  selected  for  this  study  is  a  Silicon  Inten¬ 
sified  Target  (SIT)  vidicon  detector  purchased  from  EGfcG  Princeton  Applied  Research  Corpor¬ 
ation1  1,-1 5.  The  spectral  radiance  from  approximately  200  channels  (of  the  available  S00) 
of  wavelength  in  the  range  0.60-0.80  ym  were  measured  since  only  these  channels  registered 
significant  signals. 

Once  the  spectral  radiance  data  has  been  obtained,  the  data  processor  uses  a  curve  fit¬ 
ting  routine  to  calculate  the  temperature.  For  the  purposes  of  curve  fitting,  the  spectral 
emissivity  is  assumed  to  follow  each  of  three  relations  in  turn: 


CA  *  e.»  +  €n  A 

e,  *  e„,  +  £  A  +  e  A2 
A  12  22 


(6) 

(7) 

(8) 


where  eij  are  constant  coefficients  for  the  particular  experimental  conditions.  Equation  8 
is  employed  by  Svet  for  his  multi-wavelength  technique9. 


The  curve  fitting  routine  is  a  two  step  process.  Initially,  a  value  for  the  temperature 
is  assumed  to  determine  the  emissivity  coefficients  of  the  linear  spectral  emissivity  func¬ 
tion  (Equation  6)  by  a  least  squares  method1 *.  By  repeating  this  process,  the  routine 
searches  for  a  temperature  which  will  yield  a  minimum  £lt  term.  The  temperature  that  is 
found  is  then  used  to  calculate  a  mean  emissivity  value: 


This  provides  approximations  of  temperature  and  emissivity  which  are  used  as  the  starting 
point  for  the  second  step  of  the  process. 

Next,  the  routine  uses  a  least  squares  fit  with  a  gradient  expansion  search  algorithm17 
to  find  the  emissivity  coefficients  and  the  temperature  using  each  spectral  emissivity  func¬ 
tion  above  in  turn.  The  method  involves  expanding  the  fitting  function  in  a  Taylor's  expan¬ 
sion  with  the  adjustable  parameters  of  the  spectral  emissivity  function  and  temperature, 
then  uses  a  least  squares  method  to  find  the  optimum  values  of  these  parameters.  These 
optimum  values  are  chosen  by  minimizing  the  reduced  chi-squared  value  for  the  fit  so  that 
the  relative  decrease  in  the  chi-squared  value  for  small  changes  in  the  parameters  is  re¬ 
duced  below  a  threshold  (0.01%). 


Computer  simulations 

To  demonstrate  the  improved  multi-wavelength  technique,  a  series  of  computer  simulations 
were  deviled.  Theoretical  spectral  radiance  data  was  generated  using  the  Planck  radiation 
law  (Equation  1)  for  a  given  schedule  of  temperature,  emissivity  function,  wavelength  range, 
and  number  of  data  pairs  (see  Table  1).  A  temperature  of  1800  K  was  chosen  for  most  sched¬ 
ules  since  this  is  the  approximate  melting  temperature  of  steel.  To  simulate  experimental 
noise,  the  generated  Planck  law  values  were  multiplied  by  a  rundom  number  series  of  uniform 
distribution  around  unity.  This  is  not  necessarily  the  form  of  noise  expected  in  practice, 
but  for  simplicity  the  noise  was  assumed  proportional  to  the  spectral  radiance.  The  rms 
value  of  the  difference  between  the  random  number  series  and  unity  is  used  to  describe  the 
three  simulated  experimental  noise  levels:  0.001,  2.851,  and  5.70%. 


Table  1.  Schedules  of  Parameters  for  the  Computer  Simulations 


EaiTT-nq 

T  (K) 

£,\ 

Data  Pairs 

a  Ranqe  (ym) 

faui 

170 

0.6  -  0.8 

2 

1273 

1.0 

200 

0.6  -  0.8 

3 

1800 

1.6  -  1 . 0A 

200 

0.6  -  0.8 

4 

1800 

1.0 

100 

0.6  -  0.8 

5 

1800 

1.0 

200 

0.5  -  1.5 

These 

generated  data  sets 

were  used  to  simulate 

measured  data  for  the  curve 

fitting  rou- 

tines. 

For  each  noise  level 

,  a  mean  and  standard 

deviation  for  the  adjustable 

parameters 

were  calculated  by  repeating  the  procedure  with  20  different  sets  of  data  for  each  schedule. 
The  mean  and  standard  deviation  of  the  temperature  parameter,  which  give  an  indication  of 
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the  errors  and  scatter  in  the  temperature  measurements,  are  shown  relative  to  the  actual 
temperature  in  Figures  2-5.  The  results  for  two  noise  levels  of  schedule  3  are  given  in 
Table  2. 
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Figure  2.  Computer  simulation  results  for 
the  constant  emissivity  assumption  (Equation 
9) .  The  mean  temperature  measurement  error 
relative  to  the  actual  temperature  is  above 
(note  the  change  of  scale)  with  the  relative 
standard  deviation  of  the  measurements  below. 


0.0  2.0  4.0  6 .  C 
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Figure  3.  Computer  simulation  results  for 
the  two  term  polynomial  emissivity  assump¬ 
tion  (Equation  6).  The  mean  temperature 
measurement  error  relative  to  the  actual 
temperature  is  above  with  the  relative  stan 
dard  deviation  of  the  measurements  below. 


Table  2.  Mean  Parameter  Results  of  the  Schedule  3  ( e*  ■  1.6  -  1.01)  Simulation  (Standard 
.  . Deviations  in  Parentheses) 


Rms  £  ^  Fitting 

Noise  (%)  Equation  T  (K)  e,  e,  e. 


Equation 

T  (K) 

1531 

c, 

oTTl 

Ll 

6 

1800 

1.60 

-1.00 

7 

1831 

0.94 

0.06 

-0.48 

8 

1752 

1.30 

-1.56 

9 

1864  (86) 

0.72  (0.40) 

6 

1834  (61) 

1.37  (0.98) 

-0.83  (0.93) 

7 

1855  (53) 

0.82  (0.39) 

1.12  (0.39) 

-0.47  (0.62) 

8 

1784  (98) 

0.94  (1.33) 

-1.31  (0.57) 

As  can  be  seen,  the  introduction  of  as  much  as  5.7%  rms  noise  causes  less  than  a  0.5% 
relative  mean  temperature  error,  with  less  than  a  5%  relative  standard  deviation,  for  most 
schedules  and  spectral  emissivity  assumptions.  Table  2  shows  that  the  errors  and  scatter 
tend  to  be  much  larger  for  the  emissivity  parameters  than  for  the  temperature  determination, 
which  is  typical  of  all  of  the  simulations.  This  behavior  is  as  expected  from  the  analysis 
of  Equations  4  and  5.  Also  consistent  with  this  analysis  is  the  reduction  in  most  errors 
and  scatter  for  the  lower  temperature  simulation  (Schedule  2)  due  to  the  increase  in  Y  value 
(Equation  3)  compared  to  the  higher  temperature  (Schedule  1) .  Reducing  the  number  of  data 
pairs  from  200  (Schedule  1)  to  100  (Schedule  4)  generally  caused  an  increase  in  the  mean 
error  and  standard  deviation,  as  is  expected  by  curve  fitting  statistics. 

In  the  wavelength  regions  and  temperatures  examined,  the  Planck  radiation  law  appears 
roughly  exponential  in  nature,  as  can  be  seen  in  Figure  1.  Thus,  especially  with  the  expo— 
nential  fitting  function  (Equation  8) ,  there  is  a  possible  ambiguity  between  the  effects  of 
a  change  in  the  temperature  or  emissivity  parameters  that  cannot  be  resolved  by  the  curve 

routines  and  results  in  a  correlation  between  these  parameters.  As  is  expected  from 
this,  the  mean  errors  and  standard  deviations  are  usually  greatest  for  the  exponential  emis- 
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0.0  2.0  4.0  6.0 

RMS  NOISE  053 

figure  4.  Computer  simulation  results  for 
the  three  term  polynomial  emissivity  assump¬ 
tion  (Equation  7) .  The  mean  temperature 
measurement  error  relative  to  the  actual 
temperature  is  above  with  the  relative  stan¬ 
dard  deviation  of  the  measurements  below. 
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Figure  S.  Computer  simulation  results  for 
the  two  term  exponential  emissivity  assump¬ 
tion  (Equation  8) .  The  mean  temperature 
measurement  error  relative  to  the  actual 
temperature  is  above  with  the  relative  stan¬ 
dard  deviation  of  the  measurements  below. 


sivity  function.  Extending  the  wavelength  range  to  0.5-1. 5  um  (Schedule  5)  greatly  reduces 
these  values,  especially  for  the  exponential  function,  since  the  Planck  radiation  law  curve 
does  not  approximate  a  linear  exponential  function  in  the  extended  range. 

The  generated  spectral  emissivity  values  for  Schedule  3  drop  linearly  from  1.0  to  0.8  as 
the  wavelength  increases  from  0.6  to  0.8  um.  The  constant  and  exponential  emissivity  fit¬ 
ting  functions  are  expected  to  have  difficulty  fitting  this  behavior  and  thus  cause  errors 
in  the  temperature  measurement.  However,  the  two  and  three  term  polynomial  emissivity  fit¬ 
ting  functions  would  be  expected  to  make  more  accurate  temperature  measurements.  In  fact 
aside  from  the  two  term  polynomial  emissivity  with  no  noise,  the  largest  mean  temperature 
measurement  errors  occur  with  Schedule  3.  Correlation  effects  between  the  fitting  para¬ 
meters  are  again  the  cause.  A  clear  example  of  this  is  seen  in  Table  2  comparing  the  mea¬ 
sured  polynomial  emissivity  parameters  for  data  with  no  noise.  In  the  absence  of  correla¬ 
tion  between  the  parameters,  the  results  for  the  two  and  three  term  polynomial  functions 


should  be  identical  (with  e. 


0  from  Equation  6) . 


The  effect  of  a  temperature  gradient  within  the  source  area  was  also  simulated.  For  sim¬ 
plicity,  the  source  is  assumed  to  be  divided  in  equal  parts  symmetrically  to  an  average  tem¬ 
perature  so  that  half  of  the  source  is  at  a  lower  temperature  and  half  of  the  source  is  at  a 
higher  temperature.  This  is  a  condition  which  is  more  severe  than  would  actually  exist. 
Three  gradients  are  compared  using  the  parameters  of  Schedule  1  of  Table  Is  0.0%  (1800  K) , 
1.0%  (1791  and  1809  K) ,  and  10%  (1710  and  1890  K) .  The  results  of  these  simulations  are 
shown  in  Table  3.  Even  with  this  extreme  approximation  for  the  gradient,  the  measurement 
errors  relative  to  the  average  temperature  (2. 4-3. 8%  for  the  10%  gradient)  are  less  than  the 
gradient  itself.  In  addition,  it  can  be  seen  that  even  with  thermal  gradients  on  the  order 
of  1%  in  the  measured  area,  the  fitting  routine  will  find  the  correct  average  temperature. 

Experimental  demonstration 

To  demonstrate  the  technique  experimentally,  the  temperature  of  a  resistively  heated 
platinum  strip  was  measured  using  the  MIRTMA  system.  The  spectral  emissivity  behavior  of 
this  source  was  unknown.  The  temperature  of  the  source  area  on  the  strip  was  monitored  with 
a  thermocouple  for  comparison  with  the  MIRTMA  results.  Measurements  were  made  at  four  tern- 
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peratures  within  1270-1725  K  and  the  reaulta  are  shown  in  Figure  6.  More  detailed  descrip¬ 
tions  of  these  measurements  and  the  MIRTMA  are  being  prepared  for  publication  and  are  also 
included  in  Reference  14. 
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Table  3.  Temperature  Results  (in  X)  of 
the  Thermal  Gradient  Simulations 
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Figure  6.  Comparison  of  the  MIRTMA  measure¬ 
ments  to  the  thermocouple  (T/C)  measurements 
of  the  platinum  strip  source  (perfect  agree¬ 
ment  and  5«  error  bands  are  denoted  by  solid 
and  dotted  lines  respectively) . 


The  MIRTMA  temperature  measurements  are  generally  within  5%  of  the  thermocouple  measure¬ 
ments,  and  are  within  0.5%  for  the  three  term  polynomial  emissivity  fitting  function.  The 
inability  of  the  constant  function  to  adequately  fit  the  actual  spectral  emissivity  behavior 
accounts  for  the  errors  with  that  function,  and  correlation  effects  are  the  cause  of  the 
large  errors  with  the  exponential  function. 

Summary 

This  paper  describes  an  improved  method  of  multi-wavelength  pyrometry  which  uses  curve 
fitting  techniques.  This  method  fits  multiple  spectral  radiance  measurements  to  the  Planck 
radiation  law  relation  to  determine  temperature  and  spectral  emissivity  simultaneously.  It 
assumes  only  that  a  smooth  function  exists  between  spectral  emissivity  and  wavelength  over 
the  range  of  the  observed  wavelength  band.  Relative  errors  in  the  temperature  measurements 
are  generally  an  order  of  magnitude  less  than  random  relative  measurement  errors  of  spectral 
radiance  and  than  relative  errors  in  the  calculated  spectral  emissivity  values. 

Computer  simulations  show  that  the  accuracy  and  precision  of  the  technique  increases  as 
the  wavelength  range  is  extended  toward  the  peak  in  spectral  radiance  and  as  the  number  of 
measured  data  pairs  increases.  Simulations  also  show  that  thermal  gradients  of  approximate¬ 
ly  1%  of  the  average  temperature  within  the  source  area  cause  negligible  errors  in  the  tem¬ 
perature  measurement.  The  MIRTMA  system  has  experimentally  confirmed  that  this  technique 
can  be  used  to  measure  temperatures  of  a  platinum  strip  within  1274-1724  X.  These  thermal 
measurements  agreed  with  thermocouple  measurements  to  within  5%  for  all  assumed  forms  of  the 
emissivity  fitting  function  and  within  0.5%  for  the  three  term  polynomial  emissivity  func¬ 
tion. 


Errors  may  be  caused  in  the  temperature  measurements  if  the  assumed  functional  form  of 
spectral  emissivity  with  wavelength  has  too  few  terms  to  fit,  or  is  otherwise  inappropriate 
for,  the  actual  behavior.  This  problem  is  exacerbated  by  correlation  effects  between  the 
temperature  and  emissivity  parameters  during  curve  fitting  with  the  current  procedures. 
Further  work  is  being  conducted  to  improve  the  fitting  procedure  and  minimize  correlation 
effects  between  the  adjustable  parameters.  Reducing  correlation  effects  should  increase  the 
temperature  measurement  accuracy  and  shorten  the  computation  time.  Studies  are  also  being 
conducted  to  determine  the  proper  emissivity  function  and  its  limits  of  validity 
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A  Method  of  Filming 
Metal  Transfer  in 
Welding  Arcs 

Standard  optical  components 
and  a  small  He-Ne  laser  are  used  to 
observe  metal  transfer  in  welding  arcs 


Photographic  apparatus  for  high  speed  timing  of  metal  transfer  in 
GMAW  and  SMAW  arcs.  Red  light  is  the  cotimated  laser  beam  used 


to  Ruminate  the  arc  area 


BY  C  D.  AI1EMAND,  R.  SCHOTOER,  D.  E.  RES,  AND  T.  W.  EAGAR 


During  studies  of  metal  transfer  in  gas  metal  arc  or  shielded 
metal  arc  welding,  it  is  often  useful  to  perform  high  speed 
cinematography  of  the  transfer  process;  however,  the  high 
light  emission  from  the  arc  creates  a  number  of  practical 
problems.  The  conventional  solution  is  to  flood  the  arc  area 
with  a  light  source  of  even  greater  intensity  than  the  arc.  This 
very  intense  light  is  then  photographed  through  a  very 
strong  neutral  density  filter.  While  this  approach  provides  a 
very  good  image  of  the  metal  transfer,  the  external  light 
source  is  relatively  expensive,  and  the  severe  filtering  of  the 
arc  light  can  result  in  loss  of  structural  information  about  the 
arc  plasma. 

This  article  describes  an  optical  layout  that  has  been  used 
to  take  moving  pictures  of  a  filler  metal  wire  and  droplets  in  a 
gas  metal  arc.  More  generally,  the  apparatus  is  applicable  to 
observation  of  an  opaque  object  surrounded  by  a  bright  and 
transparent  light  source  of  such  high  intensity  that  it  blinds 
the  detector  and  prohibits  direct  observation. 

This  article  describes  an  apparatus  that  produces  a  shad¬ 
owgraph  of  the  electrode  and  the  detaching  metal  drops. 
The  system  can  also  be  used  with  interchangeable  filters  that 
permit  a  portion  of  the  arc  light  to  be  transmitted,  thus 
providing  additional  information  about  the  arc  plasma. 

The  Shadowgraph 

The  0.8  mm  (0.030  in.)  diameter  output  beam  of  a  2  mW 
He-Ne  laser  is  expanded  and  spatially  filtered  into  a  20  mm 
(0  79  in.)  diameter  collimated  beam  and  used  to  illuminate 
the  object  represented  in  Fig.  1  by  an  opaque  sphere.  Since 
the  plasma  surrounding  the  object  is  transparent  at  633  nm, 
the  collimated  laser  light  passes  through  the  plasma  and  is 
concentrated  by  lens  1  (placed  at  a  distance  of  twice  its  focal 
length  from  the  sphere)  into  the  hole  of  a  diaphragm  placed 
at  the  focus  of  the  lens.  A  narrow  band  interference  filter 
which  transmits  at  the  laser  wavelength  is  placed  behind  the 
diaphragm.  All  the  laser  light  passes  through  this  filter,  except 
the  part  blocked  by  the  opaque  sphere  and  a  small  fraction 
which  is  absorbed  by  the  filter.  This  simple  apparatus  may 
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overfills  the  camera  lens  An  opaque  drop  is  shown  in  the  ob/ect 
area  and  its  shadow  is  traced  to  the  film.  The  image  position  and 
height  are  indicated  by  the  intersection  of  the  image  plane  and  field 
ray  The  chief  ray  through  the  center  of  each  imaging  lens  is 
represented  by  a  dashed  line 
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Fig  .  2-  16  mm  movie  shadowgraph  obtained  with  the  shadowgraph 
of  Figure  1.  to  which  a  Field  lens  has  been  added.  The  filer  metal 
wire  and  a  faKng  droplet  are  dearly  visible,  but  the  arc  is  not 
recorded.  The  interference  rings  are  due  partly  to  the  absence  of  a 
spatial  Filter  in  the  beam  expander  and  partly  to  the  presence  of  a 
narrow  band  interference  Fiter 


thus  be  used  to  form  a  shadowgraph  of  the  sphere.  The 
image  is  produced  at  one  focal  length  of  lens  q  behind  the 
diaphragm  and  is  relayed  to  the  high  speed  camera  film 
plane  by  lens  2  of  the  camera. 

To  understand  how  the  light  from  the  arc  is  excluded  from 
the  shadowgraph,  one  must  remember  that  light  from  the 
arc  is  not  collimated  but  is  emitted  in  all  directions.  The  arc 
light  is  not  focused  on  the  diaphragm,  and  most  of  it  is 
rejected  at  this  point.  The  small  portion  of  arc  light  that 
passes  through  the  diaphragm  (which  acts  as  a  spatial  filter)  is 
rejected  by  the  narrow  band  interference  filter.  The  arc  tight 
is  nearly  completely  blocked  by  both  the  spatial  and  the 
wavelength  filter,  and  only  the  laser  light  reaches  the  film  to 
produce  a  red  background  with  a  sharp  shadow  of  the 
object,  as  seen  in  Fig.  2. 

The  resolution  of  the  image  is  limited  by  the  diaphragm 
and  by  the  quality  of  the  imaging  optics.  Use  of  good  quality 
35  mm  camera  lenses  and  reasonably  large  diaphragm 
openings  [e.g..  1.5  to  3  mm)  have  produced  acceptable 
resolution. 

There  are  a  number  of  considerations  which  make  this 
system  more  practical: 

1)  The  object  on  the  first  image  plane  is  reversed  such 
that  the  camera  records  an  upside-down  picture.  This  is 
reversed  by  adding  another  lens  (lens  2,  Fig.  3)  that  produces 
a  secondary,  upright  image.  This  lens  can  also  be  used  as  a 
relay  lens  to  increase  the  distance  from  the  object  to  the 
camera 

2)  The  camera  lens  diameter  (lens  4,  Fig.  3)  is  usually  much 
smaller  than  the  beam  diameter  at  the  level  of  the  camera 
lens.  Therefore,  much  of  the  field  is  lost;  however,  this  may 
be  restored  to  full  object  size  by  a  field  lens  (lens  3)  placed 
on  the  second  image  plane. 

3)  The  camera  lens  has  a  diaphragm  that  is  normally  used 
to  adjust  the  film  exposure.  In  the  directed  light  system 
described,  light  from  a  point  in  the  object  area  reaches  a 
well-defined  and  small  part  of  the  camera  lens.  Therefore, 
the  camera  lens  diaphragm  does  not  act  as  a  normal 
aperture  stop,  but  rather  as  a  field  stop.  Since  this  is  not 
desirable,  the  camera  lens  diaphragm  is  kept  wide  open  and 
the  exposure  is  adjusted  by  means  of  neutral  density  filters. 
The  exposure  was  measured  roughly  by  means  of  a  35  mm 
still  camera  before  test  runs  with  the  high  speed  camera 
were  made.  The  aperture  of  the  35  mm  camera  was  opened 
wider  than  the  beam,  and  the  focal  length  of  its  lens  was 
taken  into  account  to  estimate  the  correct  exposure  for  the 
high  speed  movie. 


4)  When  studies  are  performed  in  a  glove  box,  the  laser 
beam  can  be  introduced  into  and  led  out  of  the  glove  box 
through  glass  windows,  also  shown  in  Fig.  3.  These  windows 
should  be  polished  plane  parallel  glass  tilted  about  5  degrees 


Table  1 — Equipment  Specifications 

laser: 

Spectra  Physics,  2  mW  He-Ne 
Model  145-02 

Expander  with  spatial 

Spindler  &  Hoyer  258  Model 

filter: 

038655 

lens  1: 

f  150  mm.  6  elements, 
self -designed*1 

Neutral  density  filter: 

ND  0.5.  OPCO  Labs. 

Leominster,  Mass. 

Erector  lens: 

Mamiya/Sekor  1:2/50 

Diaphragm: 

NRC  Model  101.5  opened  1.5 
mm 

Interference  filters: 

633  nm,  Andover  Corp., 
Lawrence,  Mass. 

Field  lens: 

f  200  mm  bi-converse  simple 
lens,  40  mm  diameter 

Camera  lens: 

Schneider  Tele-Xenar  1:2  8/75 

Spacer  ring: 

10  mm 

Camera: 

Redlake  Hycam,  model  41 

Optical  benches  and 
accessories: 

Klinger  Scientific 

u>Th*  I cm  may  be  repieced  by  I  commercial  t  1  SO  mm  photofrapNc  or  reproductive 
lem  with  «  deer  aperture  o(  25  mm  or  more 
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Fig  4-Gas  metal  arc  wekfng  of  titanium:  Top  Left-Cathode  plasma  jet  forms;  Top  Right-jet  repels  metal  drop;  Bottom  Left-drop 
detaches,  jet  breaks;  Bottom  Right -jet  broken.  Replacement  of  the  narrow  band  filter  'with  the  neutral  density  filer  provides  more 
information  on  the  behavior  of  the  plasma  as  we I  as  the  droplet  detachment 


from  vertical  to  prevent  a  second  surface  reflection  from 
following  the  main  optical  path. 

Addition  of  Arc  Light  Illumination 

The  shadowgraph  apparatus  produces  a  black  image  on  a 
red  background,  without  information  on  the  droplet  surface 
or  even  on  the  arc  itself.  Some  of  this  information  can  be 
retrieved  if  one  removes  the  633  nm  laser  filter,  thus  letting 
some  of  the  arc  light  through  the  system.  The  balance 
between  the  laser  and  the  arc  light  is  reached  by  adjusting 
the  diameter  of  the  diaphragm  hole  and  by  adding  a  neutral 
density  filter  where  the  narrow  band  filter  had  been  located. 


This  produces  an  image  of  the  electrode  and  droplets  which 
is  back-*uminated  by  the  red  laser  light  and  front-illuminated 
by  the  blue  arc  light  -  Fig.  4.  The  arc  itself  is  dearly  visible  and 
can  be  aBowed  to  overexpose  as  a  white  image  on  the  film 
The  instrument  specifications  of  this  system  are  given  in 
Table  1. 

Conclusion 

A  method  has  been  described  which  permits  viewing  of 
metal  transfer  in  welding  arcs  using  relatively  inexpensive 
optical  components  and  a  small  size  He-Ne  laser.  In  many 
cases,  the  optical  information  transferred  by  this  method  is 
superior  to  that  obtained  through  use  of  intense  external 
light  sources  and  very  dense  neutral  filters 
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WEAR  OF  CONTACT  TIPS  IN  GAS-METAL  ARC  WELDING  OF  TITANIUM 

by 

Jerome  Lee  LaPointe 


Submitted  to  the  Department  of  Mechanical  Engineering 
on  June  4,  1984  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

Bachelor  of  Science 


ABSTRACT 

Wear  in  copper  contact  tips  occurs  much  more  rapidly 
in  gas-metal  arc  (GMA)  welding  when  titanium  wire  is  used 
than  when  steel  wire  is  used.  This  rapid  wear  necessitates 
frequent  changing  of  contact  tips,  which  is  expensive.  Ex¬ 
periments  were  run  to  understand  the  wear  mechanisms  and 
factors  controlling  wear  in  the  contact  tip.  The  results  of 
the  experiments  reveal  that  titanium's  low  thermal  conduc¬ 
tivity  causes  melting  of  the  wire  at  the  contact  tip-wire 
interface.  Melted  titanium  can  build  up  at  the  contact 
point  and  can  lead  to  jamming  of  the  wire  in  the  contact  tip 
or  it  can  "freeze"  to  the  moving  wire  and  abrade  the  softer 
contact  tip  surface.  Jamming  of  the  wire  in  the  contact  tip 
did  not  occur  in  this  experiment. 


Thesis 
Supervisor : 
Title  s 


Thomas  W.  Eagar 

Associate  Professor  of  Materials  Engineering 
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Appendix  H 

(presented  at  AIME  Annual  Meeting,  New  York,  NY,  February  1985) 

CONTACT  TIP  WEAR  IN  GAS-METAL  ARC  WELDING  OF  TITANIUM 

by 

K.  T.  Ulrich,  J.  L.  LaPointe,  G.  B.  Hunter,  and  T.  W.  Eagar 

ABSTRACT 

When  titanium  wire  is  used  in  gas-metal  arc  (GMA)  welding,  copper  contact 
tips  wear  more  rapidly  than  when  steel  wire  is  used.  This  wear  is  accompanied 
by  arcing  between  the  base  plate  and  the  contact  tip  instead  of  between  the 
plate  and  the  filler  wire,  a  phenomenon  called  "burnback."  Burnback  requires 
a  complete  shutdown  of  the  welding  operation  and  is  therefore  an  expensive 
problem.  Experiments  were  run  to  understand  the  mechanism  of  contact  tip 
failure.  The  results  of  these  experiments  reveal  that  the  low  thermal 
conductivity  of  titanium  causes  melting  to  occur  at  the  sliding  electrical 
contact  junction  between  the  filler  wire  and  the  contact  tip.  Melting  and 
subsequent  freezing  of  titanium  at  the  contact  point  can  either  build  up  and 
cause  jamming  of  the  filler  wire  in  the  contact  tip  or  can  cause  adhesion  of 
a  chunk  of  titanium  to  the  filler  wire  that  then  abrades  the  copper  contact 
surface.  The  mechanism  leading  to  the  jamming  of  wire  in  the  contact  tip 
led  to  burnback  in  several  of  the  experimental  cases. 
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Appendix  I 
(in  press) 

Influence  of  Arc  Pressure 
on  Weld  Pool  Geometry 

A  new  model  of  a  compound  vortex  is  proposed  as  a 
possible  mechanism  to  explain  the  deep  surface  depression 
encountered  at  currents  over  300  amperes 


ABSTRACT.  At  currents  over  300 
amperes,  the  surface  of  the  wek)  pool 
becomes  marked^  depressed  and  the 
assumption  of  a  flat  surface  is  no  longer 
valid.  In  order  to  predict  the  weld  pool 
geometry,  the  shape  of  the  surface 
depression  under  the  action  of  the  arc 
pressure  has  been  calculated. 

At  currents  of  300  amperes,  it  is  found 
that  the  arc  pressure  cannot  account  for 
either  the  depth  or  the  shape  of  the 
experimental  observed  surface  depres¬ 
sion.  Instead,  a  new  model  of  a  com¬ 
pound  vortex  is  proposed  as  a  possible 
mechanism  to  explain  the  deep  surface 
depression  in  this  current  range. 

Introduction 

A  number  of  investigators  have  stud¬ 
ied  the  magnitude  of  arc  pressure  in  gas 
tungsten  arcs  (Refs.  1-6).  Some  attempts 
have  been  made  to  explain  the  formation 
of  several  weld  defects  such  as  humped 
beads,  finger  penetration  and  undercut¬ 
ting  (Refs.  6,7,8)  based  on  the  assumption 
that  the  arc  pressure  depresses  the  sur¬ 
face  of  the  weld  pool.  An  analytical 
model  was  developed  by  Friedman  (Ref. 
9)  to  simulate  the  distortion  of  fuly  pene¬ 
trated  molten  pools  in  thin  plates  under 
the  action  of  arc  pressure  and  gravitation¬ 
al  forces. 

In  our  own  experiments,  we  have 
studied  the  influence  of  welding  currents 
on  the  depth  and  shape  of  the  surface 
depression  (Ref.  10).  Surprisingly,  it  was 
found  that  the  surface  depression  depth 
is  very  smal  (less  than  1  mm,  i.e.,  0.04  in.) 


Based  on  a  piper  presented  it  the  65th  Annual 
A  WS  Convention  held  in  Dalas,  Texts,  during 
Aprl  8-13,  1984. 

M  L  UNind  T.  W  IACAR  ire  tssodited  with 
the  Materials  Processing  Center,  Massachu¬ 
setts  Institute  of  Technology,  Cimbridge,  Mas- 
sichusetts 


BY  M.  L  UN  AND  T.  W.  EAGAR 


at  currents  up  to  240  amperes  (A),  but 
the  depth  increases  rapkfy  as  the  current 
is  increased  by  30  A-Fig.  1.  This  rapid 
change  in  surface  depression  depth  is 
dtfficult  to  explain  by  tire  increase  in  arc 
pressure  which  scales  paraboficaly  with 
weld  current  (Ref.  11).  Arc  pressure  also 
wi  not  explain  why  one  current  value 
(e.g..  260  A)  can  give  either  a  depth  of  1 
mm  (0.04  in.)  or  a  depth  of  5  mm  (0.20 
in.). 

In  order  to  evaluate' whether  arc  forces 
of  the  magnitudes  measured  previously 
(Ref.  2)  are  capable  of  explaining  the 
experimentally  observed  surface  depres¬ 
sion  depths,  an  analytical  model  was 
developed.  This  model  accounts  for  the 
balance  of  hydrostatic  potential  energy 
and  surface  energy  with  the  work  per¬ 
formed  by  the  arc  pressure  displacing  the 
Squid. 

It  wB  be  seen  that  this  model  may  be 
used  to  explain  the  experimental  obser¬ 
vations  at  low  currents  where  surface 
depression  is  shalow;  however,  the 
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CIIIEIT  r  imparts  ) 

Fig.  1  -  Vitiation  of  surface  depression  depth 
with  current,  is  observed  experimentaly.  The 
arrows  denote  increasing  and  decreasing  cur¬ 
rent  and  the  fact  that  a  hysteresis  in  surface 
depression  exists  as  current  is  changed 


resuits  cannot  explain  the  experimental 
evidence  obtained  at  high  currents 
where  the  depth  is  significant.  As  a  result, 
the  assumption  that  arc  pressure  signifi¬ 
cantly  alters  weld  pool  geometry  is  no 
longer  considered  to  be  vaid.  Instead,  a 
simpified  convection  model  is  proposed 
which  can,  in  principle,  explain  the  very 
deep  surface  depression  of  the  weld 
pool  at  high  currents. 

Analytical  Model 

The  surface  depression  will  form  a 
shape  which  minimizes  the  total  energy; 
hence,  calculus  of  variations  may  be  used 
to  calculate  this  shape  under  the  action  of 
arc  pressure  subject  to  the  constraint  that 
the  volume  of  the  weld  pool  is  constant. 
In  this  model,  K  is  assumed  that  there  is  no 
convection  in  the  molten  pool  and  that 
the  weld  is  of  the  partial  joint  penetration 
type.  In  addition,  for  ease  of  analysis, 
cytndrical  symmetry  is  assumed. 

The  energy  to  be  minimized  is  the 
surface  energy  plus  the  potential  energy 
of  the  liquid  pool.  Analytically,  this  takes 
the  form 


X*  h  +  w 

,pg-2rr(h  — w)(H  — — )* 


where  the  first  integral  on  the  right  hand 
side  of  equation  (1)  is  the  interfacial 
energy  between  the  gas  and  the  liquid 
phase  and  the  second  integral  is  the 
potential  energy  of  the  liquid  pool,  with 
respect  to  a  reference  plane  at  h  *  H 
The  geometry  of  this  system  is  shown  in 
Fig.  2. 

Equation  (1)  is  subject  to  two  con¬ 
straints.  first,  constant  liquid  volume, 
which  can  be  expressed  as: 
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Rg.  2 -Geometry  of  depressed  weld  pool 
showing  the  definitions  of  several  perimeters 
used  in  this  model 


“  o  i4 

where  Gi  ■  2*rw. 

The  second  constraint  requires  that  the 
work  performed  by  the  arc  force  be 
equal  to  the  change  of  surface  energy 
and  potential  energy  of  the  Squid  pool. 
However,  the  description  of  this  con¬ 
straint  requires  some  discussion. 

The  first  law  of  thermodynamics  tels 
us  that  the  energy  of  die  system  is 
conserved;  however,  the  work  perform¬ 
ed  by  the  arc  force  is  not  a  state  function 
but  is  path  dependent.  Consider  the 
application  of  the  arc  pressure  instanta¬ 
neously.  The  volume  of  the  liquid  wX  be 
displaced  by  an  amount  AV  and  the  work 
is  PAV,  however,  this  work  is  not  revers¬ 
ible. 

The  path  of  this  form  of  work  is  given 
by  ABC  in  fig.  3  and  the  work  is  the  area 
ABCD.  If  on  the  other  hand,  one  consid¬ 
ers  that  the  pressure  is  increased  slowly 
and  incrementaly.  the  volumetric  dis¬ 
placement  of  the  Squid  might  follow  a 
path  similar  to  AC  in  Fig.  3.  The  work 
done  in  this  case  is  given  by  the  area 
ACO  This  would  be  the  reversible  work. 
If  we  assume  that  the  volumetric  dis¬ 
placement  changes  proportionally  to  the 
increased  pressure,  the  reversible  work  is 
one-half  of  the  irreversible  work  of  path 
ABC.  Depending  on  the  path  chosen,  the 
efficiency  of  work  may  vary.  Thus, 
another  variable  n  is  introduced  to 
account  for  the  efficiency  of  work  con¬ 
version. 

Accordingly,  we  may  equate  the  work 
performed  by  the  arc  force  to  the 
change  in  surface  energy  and  potential 
energy  of  the  liquid  pool  by  the  follow¬ 
ing. 


JJ2*  “  JJjdr  +  J^«dr  (3) 

where  G2  ■  ijP»c2*tw  and  JLc^dr 
is  the  fraction  of  work  done  5y  the 
arc  force  on  the  molten  pool.  The  next 
term  Gj  is  given  by. 

where  J^bjdr  is  the  difference  between 
surface  energy  with  surface  depression 
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Rg.  3— RevenUe  and  irrevenUe  PAV  work. 
ReversUe  work  is  performed  along  path  AC 
and  irrevenUe  work  is  performed  along  path 
ABC 


and  surface  energy  without  surface 
depression  of  the  molten  pool,  and  final¬ 
ly 

/-  „  _  ,  h  +  w, 

G«  “  pg2rr(h  -  w)  (H - — )  - 


performed  by  the  arc  force  and  the  total 
system  energy  which  is  the  summation  of 
potential  energy  and  surface  energy. 

If  one  apptes  the  EuieKagrangian  cri¬ 
terion  (Ref.  12), 


■—  (F  +  AiGi  +  X2Gb)  — 

vW 


a 


dr  d(dw/dr) 


(F  +  AiGn  +  XjGJ-O 


where  Xi  and  X2  are  Lagrange  mtAipiers; 
one  has  a  function  which  must  be  mini¬ 
mized  subject  to  the  two  constraints 
noted  previously. 

After  manipulation,  the  result  of  the 
Euler-Lagrangian  criterion  is 


dt 

dr 


(l+t2)* 
«rr C 


t*)»r 


pgrw+X^^-pgw, 


H  yy2 

pg2rfh(H  -  - )  ■  pg2rr(  -y  -  wH) 

where  Jjb«dr  is  the  difference  between 
potentia  energy  with  surface  depression 
and  potential  energy  without  surface 
depression  of  the  molten  pool. 

Using  equation  (2),  I  can  be  simplified 
as: 

r*  w2  h2 

J0Pg2»r(Hh  +  j  “  2 

If  we  define 

F..2w[l+(if^]*  + 

w2  h2 

pg2«r  (Hh  +  —  — ) 

then  ) «  J^dr 

Again  using  equation  (2),  equation 
(3)  can  be  reduced  to 

JS**  "  J&2dr  + 

W2 

where  Gj  »  pg2»r  —  • 


We  now  define 

G*  ■  C2  —  G3  -  G5  ■  T|Plrc2rrw  — 
<r2»r^H-(^)2J'/i  +<r2*r  — pg2*T-j- 


where 


XU 


is  the  difference  of  work 


+X1r-(1-X2)(1+t2)-^ 

,  d e  t2  I 

(fft  +  rt- — or - r)  I 

dr  1  +  t2,J 

if  one  sets 

dw 


There  are  two  boundary  conditions 
t  *  0  at  r  ■  0 

due  to  symmetry  about  the  centerline  of 
the  weld  pool  and 
w  *  0  at  r  -  R 

since  the  liquid  must  be  in  contact  with 
the  solid  surface. 

A  computer  program  was  developed 
to  evaluate  Xi  and  X2.  After  evaluating  Xn 
and  X2  together  with  suitable  boundary 
conditions  and  constraints,  equations  (4) 
and  (5)  are  then  solved  simultaneously  by 
the  Runge-Kutta  method  (Ref.  13). 

Parametric  Analysis  of  the  Model 

In  this  model,  there  are  five  parameters 
which  may  be  related  to  the  welding 
process.  These  are: 

1.  The  density  of  the  molten  metal, 
P- 

2.  The  surface  energy  of  the  molten 
metal,  a. 

3.  The  half  width  of  the  weld  pool, 
R. 

4.  The  arc  pressure,  P„c. 

5.  The  efficiency  of  work  conversion, 

The  shape  of  the  liquid-solid  boundary 
is  not  included  in  this  model.  This  is 
reasonable  if  movement  of  the  gas-liquid 
boundary  is  small  compared  to  the 
dimension  of  the  weld  pool.  It  may  be 
argued  from  hydrostatics  that  the  shape 
of  the  solid-liquid  boundary  will  not  influ- 
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Kg  4 -Assumed  surface  temperature  and  surface  tension  distribution  of  weld  pool  vs.  reduced 
weld  pool  radius.  Type  I  surface  tension  is  believed  to  correspond  to  a  steel  without  surface  active 
elements  write  Type  I  surface  tension  is  befeved  to  represent  a  steel  with  surface  active 
elements 


ence  the  result,  provided  the  gas-liquid 
boundary  does  not  impinge  upon  the 
solid-liquid  boundary. 

Steel  and  aluminum  are  chosen  in  this 
model  with  densities  equal  to  7.84  g/cm3 
and  2.7  g/cm3,  respectively.  For  the  dis¬ 
tribution  of  surface  energy  of  the  molten 
pool,  it  is  assumed  that  the  temperature 
distribution  on  the  top  surface  of  the 
molten  pool  is  as  shown  in  Fig.  4A.  The 
corresponding  surface  tension  distribu¬ 
tions  are  shown  in  Figs.  48  and  4C.  Figure 
48  is  for  pure  iron,  while  Fig.  4C  is  for 
metal  with  high  concentrations  of  surface 
active  elements,  eg.,  sulfur  and  oxygen 
(Ref.  14). 

The  magnitude  of  these  temperature 
and  surface  tension  distributions  are  not 
appropriate  for  aluminum;  however,  the 
main  purpose  of  choosing  these  distribu¬ 
tions  is  to  investigate  the  effect  of  surface 
tension  on  surface  depression  of  a  mol¬ 
ten  pool.  Other  results  could  be  calcu¬ 
lated  for  aluminum  if  accurate  surface 
tension  values  were  available.  The  half 
width  of  the  weld  pool  is  assumed  to 
range  from  0.5  to  1  cm  (0.2  to  0.39  in.). 

The  arc  pressure  is  taken  from  litera¬ 


ture  (Ref.  2)  as  shown  in  Fig.  5,  using  the 
arc  pressure  measured  at  300  A.  With 
higher  currents,  simple  multiplication  of 
the  arc  pressure  data  from  the  300  A 
curve  is  assumed.  While  this  does  not 


RADIAL  DISTANCE  C»m5 
Fig.  5 -Arc  pressure  distribution  taken  from 
literature  (Ref.  2).  Polarity -DCEN;  elec¬ 
trode-4.0  mm  (0  16  in.)  diameter  thoriated 
tungsten  (90  deg);  shielding  gas-Ar;  arc 
length  -  5  mm  (0. 197  in.) 


represent  a  physical  reality,  it  may  be 
sufficient  to  study  the  influence  of  arc 
pressure  on  surface  depression  of  the 
weld  pool. 

Two  values,  50  and  100%,  are  used  to 
account  for  die  efficiency  of  conversion 
of  work  done  by  the  arc  force  to  the 
change  of  surface  energy  and  potential 
energy  of  the  weld  pool.  The  value  of 
100%  corresponds  to  path  ABC  of  Fg.  3, 
while  the  value  of  50%  corresponds 
approximately  to  path  AC. 

Results 

The  results  of  this  model  are  shown  in 
Figs.  6  to  11.  A  more  complete  set  of 
tabulated  results  are  listed  in  Table  1.  The 
maximum  depth  of  surface  depression  of 
the  molten  pool  varies  from  0.49  to 
38.178  mm  (0.019  to  1.504  in);  however, 
for  reasonable  pressures  exerted  by  a 
300  A  arc,  the  maximum  value  is  less  than 
1.3  mm  (0.05  in.)  for  a  1  cm  (0.39  in.) 
wide  liquid  steel  pool.  This  is  much  less 
than  the  experimentally  observed  value 
of  over  4  mm  (0.16  in.)-Fig.  1. 

Density  of  Molten  Metal 

When  the  density  of  molten  metal 
increases,  the  surface  depression  depth  is 
reduced  as  shown  in  Figs.  6  to  11.  The 
effect  of  density  becomes  more  promi¬ 
nent  when  the  top  surface  of  the  weld 
pool  becomes  deeper.  In  Fig.  6,  the 
maximum  depth  of  surface  depression 
for  steel  is  0.49  mm  (0.019  in.)  and  for 
aluminum  the  maximum  depth  of  surface 
depression  is  0.642  mm  (0.025  in.) 
Hence,  reducing  the  density  by  nearly  a 
factor  of  three  increases  the  depression 
by  only  30%.  However,  at  higher  arc 
pressures,  as  seen  in  Fig.  10,  low  density 
liquid  has  a  maximum  depth  160%  larger 
than  that  of  the  heavy  liquid. 

Surface  Tension 

The  effect  of  surface  tension  on  the 
depth  of  surface  depression  of  the  weld 
pool  is  shown  in  Figs.  6  to  8  For  the 
assumed  surface  tension  distributions 
(types  I  and  II  as  shown  in  Figs.  4B  and 
4C),  the  effect  of  surface  tension  on  the 
shape  of  surface  depression  is  not  very 
important  in  both  the  low  and  the  high 
surface  depression  ranges.  The  maximum 
surface  depression  depth  is  found  to  vary 
only  20%  for  the  high  density  liquid  and 
40%  for  the  low  density  liquid  due  to 
these  different  surface  tension  distribu¬ 
tions. 

Width  of  Weld  Pool 

The  effect  of  weld  pool  width  can  be 
seen  in  Fig.  9  or  by  comparing  Figs.  6  and 
7.  When  the  width  of  the  weld  pool 
increases,  the  surface  depression  depth 
also  increases.  As  can  be  seen  in  Fig.  9, 
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fig.  6-Half  surface  depression  curve  at  300  A  arc pressue  assuming 
revers&e  work  and  a  1.0  cm  (0.39  in.)  wide  weld  pool  The  dotted 
curve  represents  an  approximate  SqukHotd  boundary 


fig.  7-Half  surface  depression  curve  at  30Q  A  arc  pressue  assuring 
reversbie  work  and  a  2.0  cm  (0.79  in.)  wide  weld  pool.  The  dotted 
curve  represents  an  approximate  fqukHotd  boundary 


when  the  surface  depression  depth  is 
smal,  the  maximum  depth  increases  only 
80%  when  the  weld  pod  width  increases 
from  10  to  20  mm  (0.39  to  0.79  in.). 
However,  when  the  surface  depression 
depth  is  large,  the  maximum  depth 
increases  140%  when  the  weld  pool 
width  increases  from  10  to  20  mm  (0.39 
to  0.79  in.). 

Arc  Premie 

The  arc  pressure  at  300  A  causes  a 
shalow  depression  (less  than  2  mm,  i.e, 
0.078  in.).  However,  when  the  arc  pres¬ 


sure  is  very  high,  as  diown  in  Figs.  10  and 
11,  the  surface  depression  becomes  very 
large.  This  may  imply  that  at  higher  cur¬ 
rents,  the  depth  of  surface  depression  is 
dominated  by  the  strong  arc  pressure. 

It  should  be  noted  that  if  the  pressure 
dstribubon  remains  fixed  in  size,  but 
increases  paraboficaly  with  current,  the 
pressures  of  2,  4  and  6  times  the  mea¬ 
sured  pressure  at  300  A  would  corre¬ 
spond  to  currents  of  425, 600  and  735  A, 
respectively.  If,  on  the  other  hand,  the 
maximum  pressure  scales  finearly  with 
current,  as  has  been  confirmed  in  our 
laboratory  (Ref.  15),  these  pressures 


wotid  correspond  to  600, 1200  and  1800 
A,  respectively. 

Work  Conversion 

The  calculated  depth  of  surface 
depression  of  irreversfcie  work  conver¬ 
sion  is  about  two  times  larger  than  that  of 
reversbie  work  conversion,  as  shown  in 
Figs.  9  and  11,  which  is  not  surprising. 

Discussion 

One  of  the  more  obvious  criticisms  of 
this  model  is  the  lack  of  a  positive  basis 


fig.  8-Half  surface  depression  curve  at  300  A  arc  pressure  assuming 
irreverside  work.  Compare  with  fig.  6.  The  dotted  cuve  represents  an 
approximate  fqukfsofd  boundary 


Fg.  10 -Half  surface  depression  curve  at  extreme  arc  pressure.  The 
dotted  curve  represents  an  approximate  SqukHotd  boundary 


fig.  9-Maximum  surface  depression  depth  vs.  weld  pod  width  for 
ddferent  materials  and  types  of  work  conversion 


Rg-  11-Maximum  surface  depression  depth  versus  arc  pressure  for 
ddferent  materials  and  types  of  work  conversion.  The  numbers  on  the 
abscissa  represent  the  mukiptcation  factor  of  the  arc  pressure  of  the 
300  A  current  in  fig.  5 
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Table  1— Calculated  Maximum  Depth  of  Surface  Depression 

Half  Max.  depth  of 


width. 

Arc 

Surface 

Work 

surface  depression, 

Material 

mmw 

pressure*"1 

tension 

conversion 

rnir 

Steel 

5 

IX 

Type  1 

Rev. 

0.490 

Steel 

S 

IX 

Type  1 

Rev. 

0.600 

Aluminum 

S 

IX 

Type  1 

Rev. 

0.642 

Aluminum 

S 

IX 

Type  * 

Rev. 

0.903 

Steel 

10 

IX 

Type  1 

Rev. 

0.893 

Steel 

10 

IX 

Type* 

Rev. 

1.029 

Aluminum 

10 

IX 

Type  1 

Rev. 

1.511 

Aluminum 

10 

IX 

Type  8 

Rev. 

1.975 

Steel 

5 

IX 

Type  1 

Irrev. 

1.000 

Steel 

5 

IX 

Type  1 

Irrev. 

1.282 

Aluminum 

5 

IX 

Type  1 

Irrev. 

1.354 

Aluminum 

5 

IX 

Type  1 

Irrev. 

1.938 

Steel 

10 

IX 

Type  1 

Irrev. 

1.846 

Steel 

10 

IX 

Type  l 

Irrev. 

2.127 

Aluminum 

10 

IX 

Type  l 

Irrev. 

3.237 

Aluminum 

10 

IX 

Type! 

Irrev. 

4.210 

Steel 

5 

2X 

Type  1 

Rev. 

1.000 

Steel 

5 

2X 

Type  8 

Rev. 

1.279 

Aluminum 

S 

2X 

Type  1 

Rev. 

1.354 

Aluminum 

5 

2X 

Type# 

Rev. 

1.938 

Steel 

10 

2X 

Type  1 

Rev. 

1.857 

Steel 

10 

2X 

Type# 

Rev. 

2.132 

Aluminum 

10 

2X 

Type  1 

Rev. 

3.235 

Aluminum 

10 

2X 

Type# 

Rev. 

4.212 

Steel 

10 

6X 

Type  1 

Irrev. 

13.162 

Steel 

10 

6X 

Type# 

Irrev. 

14.731 

Aluminum 

10 

6X 

Type  1 

Irrev. 

31.939 

Aluminum 

10 

6X 

Type# 

Irrev. 

38.178 

Steel 

5 

2X 

Type! 

Irrev. 

2.150 

Steel 

5 

2X 

Type  II 

Irrev. 

2.712 

Aluminum 

S 

2X 

Type  1 

Irrev. 

3.069 

Aluminum 

S 

2X 

Type  II 

Irrev. 

4.650 

Steel 

10 

6X 

Type  II 

Rev. 

6.821 

Aluminum 

10 

6X 

Type  II 

Rev. 

15.572 

Aluminum 

5 

4X 

Type# 

Rev. 

4.641 

Aluminum 

10 

4X 

Type  II 

Rev. 

9.372 

Steel 

7.5 

IX 

Type  1 . 

Rev. 

0.741 

Aluminum 

7.5 

IX 

Type  1 

Rev. 

1.142 

w25  4  mm  •  1  in. 

represents  mullipication  ol  the  300  A  arc  pressure  data  given  in  Fig.  5. 


for  selection  of  the  work  conversion 
efficiency,  ij.  Although  no  direct  rationale 
can  be  given  for  choosing  either  path 
ABC  or  path  AC  of  Fig.  3,  it  is  felt  that 
path  AC  is  more  realistic.  This  is  perhaps 
easiest  to  explain  if  one  considers  not  the 
total  depression  from  an  initially  flat  metal 
surface  to  a  fully  depressed  state;  but 
rather,  one  might  consider  a  small  differ¬ 
ential  depression  from  an  already 
depressed  surface,  due  to  minor  arc  pres¬ 
sure  fluctuations.  Such  a  restoration  to  a 
new  equilibrium  must  follow  path  AC 
more  closely  than  path  ABC. 

In  addition,  by  viewing  this  system  as  a 
dynamic  equilibrium  where  small  changes 
in  the  arc  force  produce  corresponding 
changes  in  the  equilibrium  surface  shap6, 
one  can  easily  understand  why  other 
forms  of  energy  {e.g.,  thermal  energy) 
need  not  be  considered  to  contribute  to 
the  formation  of  new  surface  area.  Even 
if  thermal  energy  were  to  create  new 
surface  area,  the  system  must  finally  con¬ 
sume  such  new  area  as  it  is  only  the  arc 


force  which  balances  the  surface  tension 
forces  of  the  curved  boundary.  At  final 
equilibrium,  thermal  energy  cannot  be 
considered  to  contribute  to  the  balance 
of  the  arc  force.  Only  other  forces  such 
as  surface  tension  and  hydrostatic  forces 
can  balance  the  arc  force. 

Finally,  even  if  rather  generous 
increases  in  arc  force  are  postulated,  it  is 
seen  that  the  experimentally  observed 
surface  depressions  at  300  A  cannot 
easily  be  explained  by  the  arc  force,  no 
matter  what  the  work  conversion  effi¬ 
ciency  is  assumed  to  be. 

Experimentaly,  it  has  been  observed 
that  the  surface  depression  at  currents 
below  200  A  is  negligible.  However, 
when  current  increases  to  about  300  A,  a 
deep  surface  depression  is  found,  as 
shown  in  Fig.  12.  The  maximum  depth  of 
surface  depression  is  about  4.5  mm  (0.18 
in.)  for  stainless  steel  weld  pools.  From 
the  model  of  the  arc  pressure,  it  is  seen 
that  the  arc  pressure  cannot  explain  the 
deep  surface  depression  at  300  A  cur¬ 


7-7 

fig.  12-Surface  depression  at  300  A  of  a 
stationary  arc  weld  on  Type  304  stainless  steel: 
A -top  surface  during  welting  showing  deep 
narrow  depression;  g-cross  section  of  same 
weld.  Note  the  porosity  at  the  bottom  of  this 
finger  indicating  that  the  depression  pene¬ 
trated  to  the  very  bottom  of  the  weld 


rent,  since  the  maximum  depth  of  surface 
depression  which  was  calculated  for  steel 
is  1.28  mm  (0.050  in.).  This  is  much 
smaller  than  the  actual  depth  of  surface 
depression  which  is  experimentally  found 
to  be  4.5  mm  (0.18  in.).  Thus,  some  other 
mechanism  must  be  responsible  for  the 
deep  surface  depression  at  these  higher 
currents. 

A  simplified  vortex  model  is  proposed 
in  this  work  to  account  for  the  deep 
surface  depression  at  high  currents.  In  this 
model,  a  compound  vortex  consisting  of 
a  forced  vortex  in  the  central  part  of  the 
liquid  pool  with  a  free  vortex  surrounding 
the  forced  vortex  is  assumed.*  This  is 
shown  schematically  in  Fig.  13.  The  result 
of  the  compound  vortex  model  is  shown 
in  Fig.  14,  assuming  several  angular  veloc¬ 
ities  for  the  forced  vortex.  It  can  be  seen 
that  the  depth  of  surface  depression 
caused  by  the  compound  vortex  can  be 
dose  to  the  actual  depth  of  the  surface 
depression  at  high  currents  for  angular 
velocities  between  20  and  30  rad/s.  Such 
dreumferential  flow  can  be  seen  in  high 


•A  forced  vortex  is  a  region  where  the  angular 
velocity  wi  is  a  constant.  A  free  vortex  is  one 
that  has  constant  angular  momentum,  ie  . 
Vo  •  K/r,  where  K  is  the  vortex  constant  and  9 
is  the  azimuthal  direction  of  cylindrical  coordi¬ 
nates.  The  compound  vortex  matches  a  free 
vortex  near  the  solid  boundary  with  a  forced 
vortex  in  the  center  of  the  weld  pool. 
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%  13 -Structure  of  thecampound  vortex.  Assuming  no  transition  region  between  the  forced  and 
fret  vortex 


(0.06  in.)  rarfus 


speed  movies  of  weld  pools,  which  have 
oxide  partides  floating  on  the  surface. 

Flow,  as  that  described  above,  may 
come  from  the  toroidal  flow  on  the  plane 
paralel  to  the  arc  axis.  MMere  (Ref.  16) 
found  that  in  electricaly  driven  toroidal 
flow,  a  circumferential  rotational  motion 
on  the  plane  normal  to  the  arc  axis  is 
produced  and  increases  with  the  intensity 
of  electrovortex  (toroidal  flow)  on  die 
plane  paralel  to  the  arc  axis.  Therefore, 
part  of  the  kinetic  energy  of  toroidal  flow 
on  the  plane  paralel  to  the  arc  axis  is 
transferred  to  circumferential  motion, 
inducing  flow  on  the  plane  normal  to  the 
arc  axis. 

ft  wil  be  appreciated  from  Mlere's 
work  that  such  circumferential  flow  wi 
not  be  strong  if  the  system  has  perfect 
cyfindricai  symmetry,  in  real  systems  with¬ 
out  perfect  symmetry,  the  axial  toroidal 
flow  vector  wi  be  displaced  from  the 


center  ine  resulting  in  a  circumferential 
velocity  component  which  causes  a  cir¬ 
cumferential  precession  of  the  toroidal 
flow  into  a  hefical  shape,  ft  should  be 
noted  that  the  maximum  circumferential 
velocities  assumed  in  Fig.  14  are  less  than 
one  third  of  the  velocities  noted  by 
Heiple  and  Roper  (Ref.  17)  using  high 
speed  cinematography;  hence,  these 
assumed  rotational  speeds  are  thought  to 
be  consistent  with  observed  motion  in 
actual  weld  pools. 

As  a  result  of  this  study,  three  current 
ranges  may  be  proposed  to  explain  sur¬ 
face  depression  of  weld  pools.  At  low 
currents  (below  about  200  A),  the  weld 
pool  depression  is  shallow  and  does  not 
significantly  influence  the  shape  of  the 
weld  pool.  At  higher  currents  (from 
about  300  to  500  A),  the  surface  depres¬ 
sion  greatly  influences  the  weld  pool 
geometry.  A  circumferential  vortex  flow 


may  be  responsMe  for  the  formation  of  . 
this  surface  depression  and  the  presence  i 
of  the  'finger'  penetration  seen  in  gas  . 
tungsten  arc  spot  welds,  GMAW  and 
SAW  beads.  At  even  higher  currents 
(above  about  500  A)*,  the  arc  pressure  is 
strong  enough  that  it  may  explain  some 
of  the  deep  surface  depression  which,  in 
turn,  influences  the  weld  pool  geometry. 
The  three  current  ranges  proposed  in  this 
model  are  consistent  with  the  results  of 
Qxhoski  (Ref.  18)  who  found  for  define* 
current  ranges  for  the  depth  of  penetra¬ 
tion  in  aluminum  welds.  In  his  study,  the 
penetration  shows  a  very  slow  increase 
with  current  in  the  100-200  A  range,  . 
f  otto  wed  by  a  steep  increase  with  current 
in  the  200-300  A  range.  The  300-400  A 
range  shows  almost  no  response  of  cur¬ 
rent  to  penetration,  but  in  the  450600  A 
range,  penetration  begns  to  increase 
moderately  with  current. 

Qxhoski's  100200  A  range  corre¬ 
sponds  to  the  regime  where  surface 
depression  is  controled  by  arc  force  and 
is  not  significant.  The  200300  A  range 
corresponds  to  the  transition  from  Me 
surface  depression  to  formation  of  a  fj 
vortex,  whtte  the  300-400  A  range  corre¬ 
sponds  to  bottoming  out  of  the  vortex 
depression  at  the  solid  boundary.  In  such 
a  case,  the  current  from  the  arc,  or  the 
heat,  fottows  paths  normal  to  the  isopo¬ 
tential  fines  in  the  plasma  and,  hence,  is 
concentrated  at  the  top  edges  of  the 
Squid  vortex. 


•More  recent  work  by  Lin  and  Eager  (Ref .  IS) 
has  shown  that  the  depth  of  surfs ce  depres¬ 
sion  caused  by  the  measured  ere  pressure  *t 
600  A  sid  cannot  explain  the  deep  crater  of 
the  weld  pool  found  at  even  lower  currents. 
Figure  IS  shows  the  depth  of  the  surface 
depression  caused  by  the  arc  pressure  mea¬ 
sured  in  our  laboratory. 
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Fig.  15-Half  surface  depression  curve  at  deferent  currents  assuming  reversMe  work.  These 
surface  depression  curves  are  caused  by  the  arc  pressure  measured  In  our  laboratory 
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■ig.  16  -Schematic  representation  of  the  equh 
lotential  lines  and  the  current  density  tfstribu- 
••on  when  a  deep  crater  forms.  Note  that  the 
impotential  lines  around  a  deep  depression 
jggest  that  the  arc  current  (and  hence  the 
■lermal  energy)  is  distributed  to  the  side  wats 
•t  the  pool  and  very  Sttle  current  (i.e..  beat) 
■aches  the  bottom  of  the  depression,  as  is 
spieled  by  the  shaded  region.  Increased 
• snetration  due  to  the  vortex  formation  is  thus 
•rectively  stopped 


As  such,  direct  arc  heat  does  not  reach 
1  e  bottom  of  the  vortex,  as  shown  in 
g.  16.  Hence,  no  further  melting  occurs 
t  the  bottom  of  the  depression  and  no 
irther  penetration  results.  Above  450  A, 
hi.-  arc  pressure  becomes  significant  and. 
.igether  with  the  vortex,  deepens  the 
•-pression  as  the  current  is  increased. 
Although  the  analytical  arc  force  mod- 
I  presented  here  applies  only  to  partial 
int  penetration  welds,*  it  is  believed 
■at  the  vortex  model  is  applicable  to 
>oth  complete  and  partial  joint  penetra- 
■  in  beads.  However,  more  study  of  the 
irength  and  mechanism  of  vortex  forma- 
in  is  needed  to  confirm  this  hypothesis, 
nally,  it  should  be  noted  that  if  arc  force 
.ere  in  fact  the  dominant  cause  of  sur- 
ce  depression  in  welding,  the  depres- 
on  would  change  markedly  with 
Manges  in  liquid  density. 

When  a  vortex  is  considered  the  cause 
i  surface  depression,  there  is  no  depen- 
:  -nce  of  depression  on  density.  The  fact 
iat  finger  penetrations  in  aluminum  and 
eel  are  of  nearly  equal  depth  is  consid- 
ed  supportive  of  'he  vortex  model  In 
Idition.  the  exp«-v.merul  fhape  of  fin¬ 
er  penetration,  as  seen  in  Fig.  128,  is 
iought  to  be  more  consistent  with  a 
irtex  model  than  an  arc  force  model. 


f  or  an  analysis  oi  arc  force  on  tut  penetration 
i  lds.  see  Fnedman  (Ref  9) 


'  Conclusion 

The  depth  of  surface  depression 
caused  by  the  arc  pressure  distribution 
has  been  calculated  using  calculus  of 
variations  with  suitable  constraints.  The 
results  show  that  arc  pressure  only  influ¬ 
ences  weld  pool  geometry  at  currents  in 
excess  of  500  A.  At  intermediate  cur¬ 
rents,  the  shape  of  the  pool  may  be 
influenced  markedy  by  circumferential 
convection.  This  flow  may  be  responsible 
for  the  characteristic  finger  penetration 
of  many  weld  pods.  At  even  higher 
currents,  arc  force  may  become  impor¬ 
tant  and  may  have  a  significant  influence 
on  weld  pod  geometry. 
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Refermg  to  Hg  '  1  ■»  depth  pro 

duced  by  a  compound  .ones  may  be 
e  st mated  as  tnBows 

In  the  farced  vortee  r.gpori  (r  <,  Rg). 
the  egie  vetocity  >s  constant  and 


wth  the  onpn  referred  to  ponl  a  r.  Hg 
13 

In  the  free  vortex  regno  (r  >  Ro)  the 
angular  momentum  is  constant  and 


D, 


2*  X 


with  the  onpn  referred  to  port  Bn  fg 
13  The  total  depression  D  ■  D,  4-  Dj 


UM  of  Symbols 

T  -  temperature  on  the  surface  of  weid 
pool  (°C) 

)  -  total  system  energies*) 
a  -  surface  energyferg/cm5) 
p -density  of  liquid  metaXg/cm*) 
g  -  gravitational  accelerationfcm/s2) 
h  -  ordinate  of  soM-tquki  bound- 
aryfem) 

w- ordinate  of  liquid-gas  bound- 
ary(cm) 

r  — radial  coordinate  of  weld  pooXcm) 
H- maximum  depth  of  liquid  pooKcm) 
R-  radius  of  the  top  surface  of  liquid 
pooXcm) 

i)  -  percentage  of  work  conversion(%) 
D  -  depth  of  liquid  pool  produced  by  a 
compound  vortex(cm) 

Ro-  radius  of  forced  vortex  region! cm) 
wf- angular  velocity  of  rotational 
motion  in  forced  vortex  re- 
gion(rad/s) 
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